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ABSTRACT 

Context. The APEX Telescope Large Area Survey: The Galaxy (ATLASGAL) is an unbiased continuum survey of the inner Galactic 
disk at 870 fim. It covers ±60° in Galactic longitude and aims to find all massive clumps at various stages of high-mass star formation 
in the inner Galaxy, particularly the earliest evolutionary phases. 
(ST Aims. We aim to determine properties such as the gas kinetic temperature and dynamics of new massive cold clumps found by 

ATLASGAL. Most importantly, we derived their kinematical distances from the measured line velocities. 
' ^ ' Methods. We observed the ammonia ( J,K) = (1,1) to (3,3) inversion transitions toward 862 clumps of a flux-limited sample of submm 

1 clumps detected by ATLASGAL and extracted 13 CO (1-0) spectra from the Galactic Ring Survey (GRS). We determined distances 

, for a subsample located at the tangential points (71 sources) and for 277 clumps whose near/far distance ambiguity is resolved. 

Results. Most ATLASGAL clumps are cold with rotational temperatures from 10-30 K with a median of 17 K. They have a wide 
i-C range of NH 3 linewidths (1-7 km s -1 ) with 1.9 km s -1 as median, which by far exceeds the thermal linewidth, as well as a broad 

distribution of high column densities from 10 14 to 10 16 cm" 2 (median of 2 x 10 15 crrT 2 ) with an NH3 abundance in the range of 5 to 
30 x 10~ 8 . ATLASGAL sources are massive, > 100 M , and a fraction of clumps with a broad linewidth is in virial equilibrium. We 
found an enhancement of clumps at Galactocentric radii of 4.5 and 6 kpc. The comparison of the NH3 lines as high-density probes 
with the GRS 13 CO emission as low-density envelope tracer yields broader linewidths for I3 CO than for NH 3 . The small differences in 
derived clump velocities between NH3 (representing dense core material) and 13 CO (representing more diffuse molecular cloud gas) 
1 1 ' suggests that the cores are essentially at rest relative to the surrounding giant molecular cloud. 

Conclusions. The high detection rate (87%) confirms ammonia as an excellent probe of the molecular content of the massive, cold 

■ clumps revealed by ATLASGAL. A clear trend of increasing rotational temperatures and linewidths with evolutionary stage is seen 
1 for source samples ranging from 24 fim dark clumps to clumps with embedded HII regions. The survey provides the largest ammonia 

(yQ \ sample of high-mass star forming clumps and thus presents an important repository for the characterization of statistical properties of 

. the clumps and the selection of subsamples for detailed, high-resolution follow-up studies. 
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1. Introduction and within dynamic and complex environments and are still 

deeply embedded in their earliest phases. 

. High-mass stars are essential for the evolution of galaxies be- ft is well that high-mass star formation starts in gi- 

> . cause they energize the interstellar medium and release heavy ant molecular clouds to ame e t alJ l200l) . In analogy to low- 

• ■ elements, which are determining Galactic cooling mechanisms. mass star formationj one eX pects a cold and massive starless 

. Although the understanding of high-mass star formation is con- d of > 500 Mq consisting of mo lecular g as for the ini- 

£ . sequently of great importance, only little is known about the tial conditions of massive star/c i uster formation dBeuther et all 

early stages of the formation of massive stars in contrast to ffififo . Young massive (proto) stars evolving in dense cores 

the well-established evolutionary sequence of isolated low-mass within the clumps emit ultraviolet rad i ation that leads to the 

stars (| Andreetal. || 2000|) . This is partly because of the difficult formation of u i tracompact H II regions (UCHIIRs). Many of 

observational conditions that one encounters. High-mass stars them were found in radio contin uum surveys with the VLA 

form in clusters and are rarer according to the initial mass func- ( _ | W ood & Churchwdll fl989hh . which revealed their distri- 

tion ( |Kroupaetal.| | 2011| ), therefore they are usually found at bution and number Those regions were yery helpful in tracing 

larger distances than molecular clouds, in which low-mass stars a more evolyed stage of massive star formation in me Galaxy, 

form. Moreover, high-mass stars evolve on shorter time scales often TJCHIIRs are associated with so-called hot molecular 

cores dCesaroni et al i ri992h . which are possible precursors of 

Send offprint requests to: M. Wienen UCHIIRs, also called high-mass protostellar objects (HMPOs) 

e-mail: mwienen@mpifr-bonn . mpg . de or massive young ste llar objects (MYSO s), which have also been 

* Member of the International Max Planck Research School (IMPRS) discovere d recently dBeuther et alJl2002h . Still earlier phases are 

for Astronomy and Astrophysics at the Universities of Bonn and j 1 ujjj- • 1 j u * 

_ , 1 v 3 deeply embedded in a massive envelope and are so cold that 
Cologne. 
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most ly they cannot be probed by means of mid-infrared emission 
(e.gJSridharan et al.l 120051: iPeretto & Fulled 120101: IMotte et alJ 
120071) . Under certain circumstances, these clumps can be seen 
in absorption against a bright background provided mostly 
by emission from polycyclic aromatic hyd rocarbon molecules 
(PAHs) as infrar ed dark cloud s (IRDC s, iPeraultet all Il996t 
lEgan et al.l Il998l) . ISimon et al.l (l2006bl) and IPeretto & Fuller! 
d2010t) have presented large catalogues of these IRDCs. But to 
date , only a small fraction of them has been investigated in de- 
tail dCarev et al.lll998i l2000l iPillai et al.ll2006t IRathborne et al.l 
120061) . Previous studies revealed that IRDCs are cold (< 25 K), 
dense (> 1 5 cm' 3 ) and have hig h column densities ( ~ 10 23 - 
10 25 cm 2 , IRathborne etail 120061: ISimon et al.l l2006ah . Their 
sizes (~ 2 pc) and masses (~ 10 2 -10 4 M Q ) are similar to those 
of cluster-forming molecular clumps. This, together with the 
cold temperatures and fragmente d substructure of IRDCs, sug - 
gests that they are protoclusters (IRathborne et al.ll2006l 120051) . 
Kinematic distances to IRDC samples in the first and fourth 
Galactic quadrants have been determined us ing 13 CO (1-0) and 
CS (2-1) emission ( iJackson et al.l2 006. 2008). However, surveys 
targeted at the earl y phases of high-mass star f ormation (e.g. 
IBeuther et all 120021: IWood & Churchweiil 1 1 989al) usually only 
trace one of those stages, e.g. IRDCs, HMPOs or UCHIIRs. 
In contrast, a first attempt to obtain an unbiased sample using 
a complete dust continuu m imaging at the sc ale of the molec- 
ular complex Cygnus X ( IMotte et alJ 120071) was able to de- 
rive some more systematic results on the existence of a cold 
phase (IR-quiet massive dense cores an d Class 0-like massive 
protostars) for high- mass star formation dBontemps et alj|2010t 
ICsengeri et alj|201lh . Another survey of t he molecular comple x 
NGC 6334/NGC 6357 was conducted bv lRusseil etail (1201 Ol) . 
who identified IR-quiet massive dense cores and estimated high- 
mass protostellar lifetimes. The IR-quiet massive dense cores 
are suggested to host mass ive class protostars as see n by the 
Herschel Space Telescope dNguven Luong et"aT1l201 lah . 
Only a few hundred high-mass proto- or young stellar objects 
have been studied up to now. Hence, these and especially objects 
in still earlier phases of massive star formation need to be inves- 
tigated in more detail. To achieve significant progress in that, a 
large-area survey of the Galactic plane was conducted, which 
now provides a global view of star formation at submillime- 
ter wavelengths. ATLASGAL (The APEX Telescope Large Area 
Survey of the Galaxy at 870 yon) is the first unbia sed continuum 
surve y of the whole inner Galactic disk at 870 fim ( ISchuller et al.l 
120091) . It aims to find all massive clumps that form high-mass 
stars in the inner Galaxy by using the Large APEX Bolometer 
Camera LABOCA, which is an array with 295 bolometer p ixels 
operated at the APEX telescope (ISiringo et al.ll2007L l2008h with 
a field of view of 1 1' and a beamwidth of 19.2" FWHM at the 
wavelength of 870 fim. 

ATLASGAL, which reaches a Galactic longitude of ±60° and 
latitude of ±1.5°, is able to detect objects associated with mas- 
sive star formation at various stages and to compare them 
(ISchuller et al.l 120091) . Other Galactic plane surveys, conducted 
over a similar range, are the MSX surve y covering 8 to 21 ^m 
dPrice et al.ll2001l) . 2MASS around 2 urn dSkrutskie et al. l l2006l) . 
GLIMPSE fro m 3 to 8 urn dBeniamin et al.ll2003l) . MIPSGAL at 
24 a nd 70 um dCarev et al ll2009h and HiGAL from 70 to 500 
fjm dMolinari et alj|2010l) . While submillimeter dust continuum 
surveys are essential for identifying high-mass clumps, they lack 
information on important parameters, especially the distances to 
the newly found sources. These are needed to determine other 
properties such as their masses and luminosities. But all of this 
can be addressed by molecular line observations of the high- 



mass star forming regions. Because their density is at about the 
critical NH3 density, this molecule is appropriate for determining 
properties of the clumps without contamination from large-scale 
cloud structures with lower density. S ince the molecular g as of 
the cores is very dense (~ 10 s cirT 3 IBeuther et al.l 120021) and 
cold, many molecules such as CS and CO are partly frozen onto 
dust grains. In contrast to those, the fractional g as abundance 
of am monia remains constant in prestellar cores (Tafalla et al. 
120021) . Hence, we carried out follow-up observations of northern 
ATLASGAL sources in the lowest NH3 inversion transitions. 
NH3 is known as a reliable temperature probe of inters tellar 
clouds dTafalla et al.ll2004t IWalmslev & Ungerechts|[l983T) . The 
rotational energy levels are given by the total angular momen- 
tum, J, and its projection along the molecular axis, K. Radiative 
transitions between different /^-ladders are forbidden and the 
lowest metastable energy levels, for which J = K, are thus 
collisionally excited. The intensity ratio of their inversion tran- 
sitions therefore pro vides the rotational temperature of the gas 
dHo&Townesll 19831) . which can be used to estimate the kinetic 
temperature of the cores. This is needed for proper mass esti- 
mates from the submm data. Independently, the linewidth can be 
used to derive virial masses. Moreover, the inversion transitions 
are split into distinct hyperfine components and their ratio pro- 
vides a measure of the optical depth, knowledge of which leads 
to reliable column density and rotational temperature determina- 
tions. 

As this article was ready for subm i ssion, we became aware of 
the recent study by iDunham et alJ d201 ll) . who observed NH3 
lines of sources in the first Galactic quadrant as well, identi- 
fied by the Bolocam Galactic Plane Survey (BGPS), which mea- 
sures the 1 . 1 mm conti nuum emission . In co ntrast to our sam- 
ple, the observations of IDunham et al.l d201 ll) do not cover the 
whole northern Galactic longitude range up to 60°, but are con- 
ducted only within four different ranges in longitude, and within 
a smaller Galactic latitude range of ±0.5° compared to our NH3 
survey. Our analysis is thus complementary to the results of 
IDunham etail (1201 ll) . 

In Section [2] we describe which ATLASGAL sources were se- 
lected for our NH3 follow-up observations and provide details 
of the measurements. In Section [3] we present how we reduced 
the data and derived different clump properties from fits to the 
spectra of NH3 inversion transitions. In Section[4] different cor- 
relations of the NH3 line parameters such as the velocity distri- 
bution, linewidth, rotational temperature and column density are 
shown. Then, we compare these with submillimeter dust con- 
tinuum properties of the clumps such as the H2 column den- 
sity in Section [5] and estimate gas and virial masses. We anal- 
yse additional line pa rameters derived from 13 CO (1-0) emission 
dJackson e t al. 2006) in Section|6]to investigate clump-to-cloud 
motions etc. In Section|71 the determined gas properties are com- 
pared with those of other high-mass dust-selected star forming 
regions. A summary and conclusions of this NH3 investigation 
are given in Section[8] 

This paper analyses statistics of the NH3 data obtained for the 
northern ATLASGAL sources. We provide their near and far dis- 
tance, but do not distinguish between both, although some of the 
molecular clouds have kno wn distances from a previous study 
dRoman-Duval etai 1 120091) . In a second paper (Wienen et al. in 
prep.) ammonia lines of ATLASGAL sources in the southern 
hemisphere will be studied and distances to southern molecular 
clouds will be derived. 
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2. Observations 

2.1. Source selection 

From a preliminary ATLASGAL point source catalogue, we se- 
lected a flux-limited subsample of compact (smaller than 50") 
clumps down to about 0.4 Jy/beam peak flux den sity, whose 
positi ons were selected using the Miriad task sfind dSault et al.l 
Il99l . and observed ~ 63% of these 1361 sources. From 2008 
to 2010 we observed the ammonia emission of a total of 862 
dust clumps in a Galactic longitude range from 5° to 60° and 
latitude within ±1.5°. The 870 /urn peak flux density distribution 
of observed sources is plotted as a solid black curve in Fig. Q] 
and that of the whole ATLASGAL sample of 1361 sources as 
a dashed red curve. Their comparison reveals that the 870 /im 
peak flux distributions of both samples are similar and there- 
fore we expect the NH3 subsample to be representative of the 
ATLASGAL sample as a whole. Most NH3 observations as well 
as most ATLASGAL sources from the whole sample have a peak 
flux density of ~ 0.8 Jy/beam, the distributions above those val- 
ues can be described by a power law. 

We estimate the ATLASGAL sample to be complete above 
1 Jy/beam peak flux density, which results in 644 observed 
clumps. Of those the NH3 sample contains 406 sources (63%) 
within this flux limit. Seventy-three cl umps (11%) were ob- 
served in previous ammonia surveys dS ridharan et al.l 12001 
IChurchwell et al.lll990tlMolinari et al.lfl996HPillai et al.ll2006l) . 



0.4 - 




-1 1 2 

Iog10(peak flux density/Jy beam ) 

Fig. 1. Relative number distribution of NH3 observations as solid 
black curve and of the whole northern ATLASGAL sample as 
dashed red curve vs. the logarithm of the 870 fim peak flux den- 
sity. 



2.2. Observational setup 

Observations of the NH3 (1,1), (2,2) and (3,3) inversion transi- 
tions were made from 2008 to 2010 with the Effelsberg 100-m 
telescope. The frontend was a 1.3 cm cooled HEMT receiver, 
which is used in a frequency range between 1 8 GHz and 26 GHz 
and mounted at the primary focus. In 2008 January, two dif- 
ferent spectrometers were used, the 8192 channel autocorrelator 
(AK 90) and the Fast Fourier Transform Spectrometer (FFTS). 
In May 2008, the AK90 was replaced by the FFTS, which pos- 
sesses a greater bandwidth. The autocorrelator contains eight 
individual modules. A bandwidth of 20 MHz was chosen for 
each of them, which results in a spectral resolution of about 0.5 
km s , while the FFTS consists of two modules with a cho- 
sen bandwidth of 500 MHz each and a spectral resolution of 



~ 0.7 km s . Both spectrometers are able to measure two po- 
larizations of the three observed NH3 lines simultaneously. The 
scans observed in 2008 January with the two spectrometers were 
summed to obtain a better signal-to-noise ratio. The beamwidth 
(FWHM) at the frequencies of the (1,1) to (3,3) ammonia inver- 
sion transitions at about 24 GHz is 40". Pointed observations 
towards the dust emission peaks were conducted in frequency- 
switching mode with a frequency throw of 7.5 MHz. The median 
system temperature was about 70 K, with a zenith optical depth 
that usually varied between 0.02 and 0.09 in the beginning of the 
year in about January and was more constant, about 0.07 some 
months later, in about May and October. The total integration 
time for each source was about 5 minutes. For some clumps at 
low declination, which could only be observed at low elevations 
and consequently have high system temperatures of ~ 150 K, 
longer integration times of up to 25 minutes were chosen. 

3. Data reduction and analysis 

To reduce the NH3 spectra we used the CLASS softwar^B 
Because the ammonia lines from each source are detected in 
two polarizations, all scans that belong to one polarization of an 
inversion transition were summed and then the spectra of both 
polarizations were averaged. The frequency switching results in 
two lines, which were inverted and separated in frequency by 7.5 
MHz. This procedure was reversed by a folding, which shifts the 
two lines by 3.75 MHz and subtracts one of them from the other. 
The large frequency throw led to fluctuations of the baseline, 
which needed to be corrected. For this we subtracted a polyno- 
mial baseline of the order of 3 to 7, which produces inaccura- 
cies in the subsequent analysis of the spectra. The NHs(l,l) hy- 
perfine structure of very many ammonia sources shows different 
line profiles with various linewidths, which had to be reduced in 
the same way. Hence, we developed a strategy to set windows, 
spaces around NH3 lines, that were excluded from the baseline 
fitting: The width of the window around each hyperfine com- 
ponent of the (1,1) line was chosen according to the width of 
the main line; the order of the baseline depends on the linewidth 
as well. A polynomial of high order, up to 7, was subtracted for 
lines with narrow width, < 2 km s _1 , while a low order, mostly 3, 
was chosen for lines with a broad width, greater than 4 km s _1 . 
To test if this introduced any systematic errors, we compared 
the linewidths of simulated NH3 (1,1) lines with the linewidths 
resulting from a fit of the same spectra added on observed base- 
lines. For linewidths more narrow than 1.5 km s the model 
and the fitted observation agreed, the deviation of the two be- 
comes sligthly higher with increasing linewidths. However, the 
contribution of systematic errors for lines with broad linewidths 
is only small, ~ 6% on average. Some examples for reduced and 
calibrated spectra of observed (1,1) to (3,3) inversion transitions 
are given in Fig. [2] 

The hyperfine structure of the NHKLl) line was fitted with 
CLASS. Assuming that all hyperfine components have the same 
excitation temperature, a nonlinear least-squares fit was made to 
the observed spectra with the CLASS software, which yielded as 
independent parameters radial velocity, v'lsr, the linewidth, Av, 
at the full width at half maximum of a Gaussian profile and the 
optical depth of the main line, r m , with their errors. 
Because the hyperfine satellite lines of the (2,2) and (3,3) tran- 
sitions are mostly too weak to be detected, their optical depth 
could not be determined. A single Gaussian was fitted to the 
main line of the (2,2) and (3,3) transitions with CLASS. The 

1 available at |http://www.iram. fr/IRAMFR/GI LDASl 
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three independent fit parameters were the integrated area, A, the 
radial velocity of the main line, vlsr, and the FWHM linewidth, 
Av. The temperature of the NH3 (1,1) line was derived from the 
peak intensity of the Gaussian fit to the main line. For the (2,2) 
and (3,3) transitions the temperature were calculated with 



with the error given by the rms value. If lines were not detected, 
the noise level was used as an upper limit. We list the positions 
of the sources, the optical depths of the NH3 (1,1) lines, t(1,1), 
their LSR velocities, v(l,l), linewidths, Av(l,l), and main beam 
brightness temperatures, Tmb(^ >1)> together with their formal 
fit errors in Table [T] The LSR velocities, linewidths, and main 
beam brightness temperatures of the (2,2) and (3,3) lines to- 
gether with their formal fit errors are given in Table [2] The other 
physical parameters such as the rotational temperature (T m t), 
the kinetic temperature (7\i n ) and the ammonia column den- 
sity (jVmjj ) were derived using the standard formulation for 
NH 3 spectra dHo & Towneslll983t luiigerechts et alj|1986l) . see 
Sect. 14.31 They are given with the errors (lcr) calculated from 
Gaussian error propagation in Table [3] 

3. 1 . Line parameters derived from ammonia 

Of the 862 observed sources, 752 were detected in NH3 (1,1) 
(87%), 710 clumps (82%) in the (2,2) line and 415 sources 
(48%) in NH 3 (3,3) with a signal-to-noise (S/N) ratio > 3. The 
hyperfine structure of most (1,1) lines is clearly detected, while 
that of the (2,2) and (3,3) lines is mostly too weak to be visible. 
The sample contains sources of different intensities; examples of 
them are shown in Fig. [2] There are many strong sources such as 
G5. 62-0. 08 with an S/N ratio of 30. An example for a weak line 
is the (3,3) line in G10.99-0.08. Its signal-to-noise ratio is 3, its 
fit is shown together with that of some other sources in Fig. [2] 
The comparison of peak intensities using Tables Q] and [2] yields 
an average peak intensity of the (2,2) line of 53% of the average 
(1,1) peak line intensity and an average (3,3) main beam bright- 
ness temperature of 40% of the average (1,1) peak intensity. 
Moreover, a few sources like G10.36-0.15 and G23.27-0.26 
(see Fig. |2| exhibit two different velocity components. Those 
are probably two NH3 clumps, that lie on the same line of sight 
to the observer with different velocities, e.g. 1 1 km s _1 and 44 
km s _1 for G10.36-0.15, at different distances. The two main 
components of the inversion transitions can be separated, while 
the satellite lines partly interfere with each other. These sources 
were fitted with a model with two velocity components. 

4. Results and analysis 

We first plot different correlations of line parameters, which are 
derived from the NH3 spectra observed with the Effelsberg 100 
m telescope (FWHM beamwidth = 40") (Fig. [3] - [13). An am- 
monia (3,3) maser line together with the thermal emission is 
displayed in Fig. [14] Then, ammonia results are compared with 
the 870 yt/m flux integrated over the Effelsberg beam, extracted 
from s ubmillimeter dust co ntinuum maps of the ATLASGAL 
project dSchuller et alJl2009h (Figs.[[5]and[T6]), and with line pa- 
rameters from the 13 CO (1-0) emission detected by the Boston 
University-Five Co llege Radio Astrono my Observatory Galactic 
Ring Survey (GRS lJackson et al.ll2006l) . Fig.fTTl compares NH3 
(1,1) and 13 CO (1-0) lines. Fig. [18] shows our analysis of 13 CO 
and NH3 line-centre velocities, both molecular lines for extreme 



sources of this plot are displayed in Fig. [19] NH3 (1,1) and 13 CO 
(1-0) linewidths are compared in Fig.l20l 

4. 1 . Ammonia velocities 

The distribution of NH3 velocities of the ATLASGAL clumps 
with Galactic longitude is shown in Fig. [3] For most clumps, the 
velocities of the ammonia (1,1) line lie between -10 km s _1 and 
120 km s -1 , only few exhibit extreme velocities of between -20 
km s" 1 and -30 km s _1 and between 120 km s to 150 km s . 
They are compared w ith CO (1-0) emis sion, observed with the 
Cf A 1.2 m telescope ( iDame et"al Bool , which is shown in the 
background. Most clumps are strongly correlated with CO emis- 
sion, which probes the larger giant molecular clouds, but some 
with extreme velocities seem to be related to only weak CO 
emission. These probably consist of compact clumps, which are 
detected with the Effelsberg beamwidth of 40", but are not visi- 
ble in the larger CfA telescope beam of 8.4'. The straight line in- 
dicates the 5 kpc molecular ring that represents the most massive 
concentratio n of molecular gas an d star formation activity in the 
Milky Way (ISimon et al.ll2006bl) . Its high CO intensity results 
from integration over many molecular clouds, which are iden- 
tified by 13 CO obse rvations of the Galactic Ring Survey (GRS, 
ISimon et alj|2006bl) . For the ATLASGAL sample the NH 3 ve- 
locities are crucial to obtain kinematical distances to the clumps, 
which were calculated usin g the revised rotat ion parameters of 
the Milky Way presented bv lReid et all d2009l) . We have derived 
near and far distances (cf. Table HI and we will distinguish be- 
tween them by comparing ATLASGAL with extinction maps 
and absorption in HI lines (Wienen et al. in prep.). Sources as- 
sociated with infrared extinction peaks are most likely at the 
near distance. Hence, the addition of the ammonia measurements 
to the ATLASGAL data reveals the three-dimensional distribu- 
tion of massive star forming clouds in the first quadrant of our 
Galaxy. Some clumps are located within GRS molecular clouds, 
whose kinematic distances have already been determined by us- 
ing absor ption features in the HI 21 cm line against continuum 
emission dRoman-Duval et al I2009I These known distances are 
marked by a star in Table|4] In addition, sources located near the 
tangential points have similar near and far distances, which are 
marked by a star in Table |4] 

The distribution of Galactocentric distances can be studied in- 
dependently of near/far kinematic distance ambiguity. Fig. |4] 
shows an enhancement of the number of ATLASGAL sources 
at Galactocentric radii of 4.5 and 6 kpc, where the Scutum arm 
and t he Sagittarius arm are l ocated. This agrees with previous re- 
sults. [Ng]iyenTuon£etaI] J20TT3) used ATLASGAL and com- 
plimentary data to show that the region around Galactic longi- 
tude of 30°, where the W43 Molecular Complex is located, has a 
very high-mass concentration of sources and star formation ac- 
tivity at a Galactocent ric radius of 4.5 kpc. Th e two peaks in Fig. 
3] are also detected in lBronfman et al.l d2000l) . based on a CS(2- 
1 ) survey of IRAS point-like sources with FIR colours of UCHII 
regions (see Fig. 3 and Table 2). Moreover, the same features 
of Galactic structure were discovered already by the first H II 
region radio recomb ination line survey dReifenstein et al. |[197§ 
and more recently in lAnderson & Banial d2009l) in a sample of H 
II regions over the extent of the GRS (see Fig. 9). 

4.2. Linewidth 

The NH3 (1,1) inversion transitions exhibit widths from 0.7 
km s~ ! to 6.5 km s _1 (cf. Fig. [5] upper panel), while those of 
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Fig. 2. Reduced and calibrated spectra of observed NH3(1,1), (2,2) and (3,3) inversion transitions, the fit is shown in green for some 
sources. Towards G10.36-0.15 and G23.27-0.26 two lines are seen at different LSR velocities. 



the (2,2) lines range up to 7.5 km s _1 . Linewidths narrower than 
0.7 km s _1 cannot be resolved with the backend spectral resolu- 
tion. One of the clumps with broad linewidths is G10.47+0.03 
with Av(l,l) = 6.4 km s an d Av(2,2) - 6.7 km s " 1 , which is a 
bright, well-known hot core (ICesaroni et al I ll992h . Because the 
hyperfine structure of the (2,2) line of most clumps could not 
be detected, we determined the linewidth by a Gaussian fit. The 
NH3(1,1) and (2,2) linewidths are compared in the lower panel 
of Fig. [5] the straight dashed line in green shows equal values. In 



addition, we created a contour plot by dividing the (1,1) and (2,2) 
linewidth ranges into bins of 0.1 km s _1 and counted the num- 
ber of sources in each bin (upper panel of Fig. |5J. The widths of 
the (2,2) transition, fitted by a Gaussian, are on average broader 
than those of the (1,1) line. The fitting of the (2,2) transition by 
a Gaussian does not take optical depths effects and the underly- 
ing hyperfine structure splitting into account. For the (2,2) line 
the Gaussian fits can therefore result in broader linewidths than 
for the (1,1) transition, for which we used a hyperfine structure 
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non-thermal contribution. The thermal linewidth is given by 



40 30 20 

Galactic Longitude (deg) 



Fig. 3. LSR velocities of detected sources ar e plotted against th e 
Galactic longitude with CO (1-0) emission ( iDame et al.ll2001l) . 
which is shown in the background. The st raight line denotes the 
5 kpc molecular ring dSimon et al.ll2006bb . 



100 




galactocentric distance (kpc) 



Fig. 4. Number distribution of northern ATLASGAL sources 
with Galactocentric distance. 



fit. To investigate this difference, we used a subsample of south- 
ern ATLASGAL sources with higher S/N data (Wienen et al. 
in prep.), for which we were able to detect the hyperfine struc- 
ture of the (2,2) transition, and compared their linewidths de- 
rived from Gaussian and hyperfine structure fits. This resulted in 
the correlation Av'Gauss 

(2,2) = 1.15- Av H fs(2,2) + 0.13, which 
is a better fit to the data as shown by the solid red line in Fig. 
|5] The contribution to the line broadening from the optical depth 
of the line is 15% with a (2,2) optical depth of 1.24, an addi- 
tional broadening results from the magnetic hyperfine structure 
of the (2,2) line, which is only of about 0.1 km s _1 and there- 
fore cannot be resolved. The contour plot in the top panel of Fig. 
[5] shows that the contour lines are distributed equally around the 
red line. However, the lower panel reveals that some sources still 
lie above the red line, which indicates broader (2,2) than (1,1) 
linewidths. Because we assumed that the same volume of gas 
emits the (1,1) and (2,2) inversion lines, the same beam filling 
factor was used for all transitions in calculating the temperature. 
The difference in the (1,1) and (2,2) linewidths now shows that 
some measured NH3 lines do not exactly trace the same gas and 
using equal beam filling factors for those is therefore only an ap- 
proximation. 

The widths of the NH3 inversion lines consist of a thermal and a 



Av, 



therm — 



8^(2)^ 



'«NH ; 



(2) 



with the kinetic temperature T^ n , the Boltzmann constant k and 
the mass of the ammonia molecule twnh, ■ Assuming a mean T^ n 
of 20 K, the thermal linewidth is about 0.22 km s _I . On average 
measured NH?(1,1) linewidths are ~ 2 km s , i.e. they are dom- 
inated by non-the rmal contributions . 

In comparison to iJiiina et al] d!999l) . who measured NH3 lines 
of low-mass cores and derived a mean linewidth of about 0.74 
km s _1 , the linewidths of the ATLASGAL sources are broader 
on average. The beamwidth of the Effelsberg 100 m telescope is 
40" at frequencies of the NH3 lines, corresponding to a linear 
scale of about 0.8 pc at a typical near distance of our sources 
of ~ 4 kpc. Hence, several cores might fill the Effelsberg beam 
(see discussion of beam filling factors in Sect. 14.51 ). Assuming 
that each core possesses a narrow linewidth, they may add up to 
the observed linewidths due to a velocity dispersion of the cores 
within the telescope beam. High-resolution measurements using 
interferometers are important for investigating individual cores, 
because they might reveal possible non-thermal contributions to 
the linewidth of each of them. 

To characterize the ammonia sample, Figs.|6]and|7]illustrate the 
main beam brightness temperature and the optical depth of the 
NH3 (1,1) lines plotted against the (1,1) linewidth. The binning 
used for the (1,1) linewidth is 1 km/s, the (1,1) peak intensity 
values are divided into bins of 1 K and the (1,1) optical depth 
into bins of 1 as well. The (1,1) peak intensity ranges between 
0.2 and 6 K with a peak at 1 K, as displayed by the contour plot. 
Values of the (1,1) optical depth lie between 0.5 and 5.5 with 
an average error of 14% and a peak at 1.8. There are no clear 
trends seen in the (1,1) main beam brightness temperature and 
the optical depth. 

4.3. Rotational temperature 

The rotational temperature between the (1,1) a nd (2,2) inver- 
sion t ransition can be determined via the relation (iHo & TownesI 
119831) 



-41.5 



^(l-ex P (-T m( l,l))))) 



(3) 



with the optical depth of the (1,1) main line, r m (l, 1), and the 
main beam brightness temperatures, 7mb, of the (1,1) and (2,2) 
inversion transitions. The NH3 lines of most sources are well fit- 
ted, although the maximum intensity of some lines is slightly 
underestimated by the fit. The rotational temperature lies be- 
tween 10 K and 28 K with an average error of 9% (Table [3]) 
and a peak at ~ 16 K; it is plotted against the (1,1) linewidth in 
Fig. [8] We chose only sources with reasonably well-determined 
rotational temperatures (errors smaller than ~ 50%) for the cor- 
relation plot. There is a trend of increasing rotational tempera- 
ture with broader width of the (1,1) lines. This, together with 
the broadening of the (2,2) lines, which probe higher tempera- 
ture regions compared to the (1,1) lines, shows an increase of 
turbulence with temperature and is an indication for star forma- 
tion feedback such as stellar winds or outflows. In contrast, some 
clumps exhibit similar (2,2) and (1,1) linewidths and low rota- 
tional temperatures. Those sources could be prestellar clumps 
without any star formation yet and therefore low temperatures. 
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Av(1 ,1) (km/s) 



Av(1 ,1 ) (km/s) 




Av(1,1) (km/s) 

Fig. 5. Correlation plot of the NH 3 (1,1) and (2,2) linewidths. 
The straight green dashed line corresponds to equal widths, the 
red solid line includes the correction for the hyperfine structure 
of the (2,2) line Av Gauss (2,2) = 1.15 • Av HFS (2,2) + 0.13. The 
(1,1) and (2,2) linewidths are mostly equally distributed around 
the red solid line, as shown by the contour plot in the top panel. 
For this upper plot we counted the number of sources in each 
(1,1) and (2,2) linewidth bin of 0.1 km s . The contours give 10 
to 90% in steps of 10% of the peak source number per bin, these 
levels are used for all contour plots in this article. 



*t 



n-fn m 




Av(1 ,1) (km/s) 

Fig. 6. NH3 (1,1) main beam brightness temperature plotted 
against the (1,1) linewidth. The upper panel displays the con- 
tour plot with a binning of 1 km/s for the (1,1) linewidth and of 
1 K for the (1,1) temperature, while the lower panel shows the 
scatter plot. 



and r ro t of 22 K, it is also asso ciated with an IRAS point source 
dZoonematkermani et al.|[T990l) . 



There are also some clumps with broad (1,1) lines and low ro- 
tational temperatures that have large virial masses. The virial 
mass is a measure of the kinetic energy of a cloud (cf. Sect. 15. 21 ). 
Some examples of the most extreme sources in Fig.[8]are, e.g. the 
clump with an extremely narrow (1,1) linewidth and/or low rota- 
tional temperature, G34.37- 0.66 with Av(l,l) - 0.7 8 km s -1 and 
r rot = 10 K, a known IRDC dPeretto & Fullerll2009l) . Among the 
clumps with broad linewidth are we ll-known sources, su ch as 
the UCHIIR/hot core G10.47+0.03 dCesaroni et all 1 19921) (see 
Sect. 03} with Av(l,l) = 6.38 km s _1 and r rot = 21 K. Another 
source with a high T ro t of 26 K and Av(l,l) of 5.06 km s" 1 is 
G43 .8 0-0. 13, well-known for h arbouring OH dBraz & Epchteinl 
[1981 and H2O maser emission dLekhtll2000l) . G12 .21-0.10has 
alread y been observed with the VLA at 5 GHz by iBecker et al.l 
(11994 . who identified it as a possible UCHII region. Hence, this 
clump is probably in a more evolved phase of high-mass star for- 
mation and therefore has a high Av(l,l) of 5.73 km s _1 as well 



4.4. Source-averaged NH 3 column density and kinetic 
temperature 

To calculate the source-averaged ammonia column density, the 
optical depth and linewidth of the (1,1) inversion transition are 
needed as well as the rotational temperature, which is derived 
from equation [3] Because the optical depth and the rotational 
temperature depend only on line ratios, the resulting column 
density is a source-averaged quantity. The values of the rota- 
tional temperature, kinetic temperature, and the logarithm of the 
NH3 column density are given in Table [3] The total column 
density was derived from the column density of the (1,1) level 
(iV(l, 1)) assuming that the energy levels are populated via a 
Boltzmann distribution, defined by the rotational temperature. 
Another assumption is that only the four lowest metastable lev- 
els are populated because of the low rotational temperatures of 
the cloud clump, which is ranging from 10 to 25 K, and H2 den- 
sities of 10 s cirT 3 . Hence, the total column density is given by 
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Fig. 7. NH 3 (1,1) optical depth compared with the (1,1) 
linewidth. Contour lines of both parameters are plotted in the 
upper panel with the range of (1,1) linewidths divided into bins 
of 1 km/s and that of (1,1) optical depths into bins of 1. 



lRohlfs& Wilsonl(l200l 



3 6XP 



N tot * N(l, 1) 



with 



N(l, 1) = 2.07 ■ 10 



23.1 



r ro t(i,2) 



5 / 41.2 

1 + r xp hwi^) 



(4) 



14 
+ y exp 



99.4 

>rotd,2) 



gl v- T m t 
guAul 



rAv, 



(5) 



where gi and g u are the statistical weights of the lower and up- 
per inversion level, A u \ the Einstein A coefficient, v the inversion 
frequency in GHz and Av the linewidth in km s . Weak sources, 
whose errors of the column density are greater than 50%, were 
excluded from correlation plots (cf. Fig. [9]). 
The kinetic temperature is obtained from the rotational tempera- 
ture dTafalla et all2004l) : 



T\dn = 



7rot(l,2) 



1 - 



3t(L2) 



ln(l + l-lexp(^))' 



(6) 



where the energy differ ence between the NH 3 (1,1) and (2,2) in- 
version levels is 42 K. iTafalla et alJ (120041) ran different Monte 



Av(1 ,1) (km/s) 

Fig. 8. Dependence of the rotational temperature between the 
(1,1) and (2,2) inversion transition on the NH3(1,1) linewidth 
is shown as contour plot in the top panel and as a scatter plot in 
the lower panel . The binning of the (1,1) linewidth in the contour 
plot is 1 km/s and that of the rotational temperature is 2 K. 



Carlo models to compare their derived rotational temperatures 
with the kinetic temperatures used in the models. Equation[6]is 
derived from their fit of the kinetic and rotational temperatures. 
These authors predicted kinetic temperatures in the range of 5 K 
to 20 K, in which most of our sources fall with a peak of T^j, 
at 19 K, to better than 5%. The logarithm of the averaged NH3 
column density is compared with the (1,1) linewidth in Fig. [9] 
and with the kinetic temperature in Fig. [10] The column den- 
sity ranges between 4 x 10 14 cm -2 and 10 16 crcT 2 with an av- 
erage error of 17% (cf. Table [3) and a peak at 2 x 10 15 cm" 2 . 
The source w ith the highest colu mn density of 10 16 crrT 2 is 
G10.47+0.03 dCesaroni etal.119921) . The range of column densi- 
ties of ATLA SGAL sources compares well with values found for 
UCHIIRs bv IWood & Churchwelll (Il989bl) and high-mass pro- 
tostellar objects (H MPOs), that are considered to be in a pre - 
UCHII region phase (iMolinari et al.ll2002t iBeufher et al]l2002l) . 
It is som ewhat above that o f low-mass, dense cores from a NH3 
database dJiiina et al 119991 whose column densities lie between 
10 14 cur 2 and 3 x 10 15 cm 2 . 
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log10(N NH /cm ) 




log10(N NH3 /cm ") 

Fig. 9. Correlation plot of the logarithm of the column density 
and the (1,1) linewidth. The contour plot is illustrated in the up- 
per panel. For the logarithm of the column density bins of 0.4 
cm -2 and for the (1,1) linewidth bins of 1 km/s are chosen. 



4.5. Beam filling factor 

The beam filling factor, 77, which gives the fraction of the beam 
filled by the observed source, is derived from equations of radia- 
tive transfer 



WM) 



CTex-:Tbg)-(l-exp(-T(l,l)))' 



(7) 



where 7mb(1, 1) is the the main beam brightness temperature 
of the (1,1) inversion transition, t(1, 1) the optical depth of the 
(1,1) main line and 7^=2.73 K. Because the H2 density of 
the ATLASGAL clumps is ~ 10 5 cnr 3 dBeuther etal.ll2002t 
iMotte et all 120031) . we assumed local thermodynamic equilib- 
rium (LTE) and accordingly used the kinetic temperature as ex- 
citation temperature in Eq. [7] We found filling factors mostly in 
the range between 0.02 and 0.6 with one higher value of ~ 0.9 
(cf. Table |4]i. Their distribution is shown as a histogram in Fig. 
[TTI Equation [7] can lead to underestimation of the filling factor 
in cases of sub-thermal excitation. To investigate the deviation 
of the excitation temperature from the kinetic temperature, we 
used the non LTE molec ular radiative transfer program RADEX 
(Ivan der Tak et al.l|2007l) . For a background temperature of 2.73 
K, an average kinetic temperature of our sample of 19 K, an 
H2 density of 10 5 cirT 3 , a column density of 10 15 cm~ 2 and a 




15 15.5 16 

log10(N NH /crrr 2 ) 

Fig. 10. Logarithm of the column density compared with the ki- 
netic temperature. There is no correlation between both. Contour 
lines of the NH3 parameters are shown in the upper plot. Bins of 
0.4 cirT 2 are used for the logarithm of the column density and 
bins of 2 K for the kinetic temperature. 



linewidth of 1 km s _1 we obtained for para-NH3 an excitation 
temperature between the lower and upper (1,1) inversion level 
of ~ 18 K. Hence, the excitation and kinetic temperature agree 
within 94%, which does not hint at sub-thermal excitation. Since 
typical ATLASGAL source sizes, obtained from Gaussian fits 
of the 870 /mi dust continuum flux, are ~ 40" (cf. Sect. 15.2b . 
low beam filling factors indicate dumpiness on smaller scales, 
as is also evident from interferometer NH3 observati ons (e.g. 
iDevine et al EoiltlRagan etd]l20Tltlolmi et al.ll20loh . 

4.6. NH 3 (3,3) line parameters 

In Fig.[T2]we plot the widths of the ammonia (3,3) line, Av(3,3), 
against the (1,1) line, Av(l,l). The straight line shows equal 
widths. A binning of 0.7 km/s is chosen for Av(l,l) and of 1 
km/s for Av(3,3) for the contour plot. While the (1,1) linewidths 
lie between 1 and 6 km s _I , the widths of the (3,3) line range 
from 1 to 10 km s _1 (cf. Table |2]i and, with a peak at 3.5 
km s , are broader than those of the (1,1) inversion transition 
for most clumps. Because the excitation of the (3,3) inversion 
transition requires a high temperature, the (3,3) emitting gas is 
likely heated by embedded already formed or forming stars that 
generate the turbulence observed in the increased linewidth. 
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Fig. 11. Histogram of beam filling factor computed according to 
equation[7] Most sources have a beam filling factor of 0.1. 



Furthermore, we investigated whether the observed NH3 (3,3) 
temperatures are consistent with the temperatures found from 
the (1,1) and (2,2) lines or whether they probe an additional, 
embedded hot core component of the clumps. We calculated the 
expected envelope (3,3) main beam brightness temperature from 
equations of radiative transfer under LTE conditions using the 
results from Sects. 1431 and 14.41 and compared it to the measured 
?mb(3,3). 

Sources towards which the NH3 (3,3) line is not detected have an 
average (1,1) linewidth of 1.8 km s _1 , which is slightly below the 
mean value of the whole ATLASGAL sample of 2 km s _1 . Their 
average rotational temperature is 15.3 K, which is also lower 
than the mean of 16.5 K. 

Assuming LTE (and thus 7ki n =r ex ), the (3,3) main beam bright- 
ness temperature is 



r MB (3, 3) ca i c = 77 • r kin (l - exp(-r(3, 3))) 



(8) 



assuming the same beam filling factor as for the (1,1) line (see 
Eq.|7]i and the optical depth t(3, 3) of the NH3(3,3) line given by 



r(3, 3) = ax 1 - exp 



hv(3,3) 
kT\ 



kin 



(9) 



where v(3, 3) is the frequency of the (3,3) inversion transition 
(23.87 GHz), k the Boltzmann constant and 7\j n the measured 
kinetic temperatures, a is defined as 



2 Vm2 1 



JV(3,3)A M/ (3,3). 



8ttv(3,3) 3 Av(1,1) 

Here, c denotes the speed of light, Av(l,l) the observed (1,1) 
linewidth, N(3, 3) the (3,3) level population number and A„/(3, 3) 
the Einstein A coefficient of the (3,3) transition. From the 
Boltzmann distribution of energy levels we determined the (3,3) 
level pop ulation number using de rived column densities Ntot (cf. 
Table I3J dRohlfs & Wilsonll2004l) 

A(3,3) = AWe, (10) 
where Q is given by 

„ 1 /123.57\ 3 /99.4\ 5 /59.12\ , 
(J - — exp H exp H exp + 1 . 

The ratio of observed to calculated (3,3) main beam bright- 
ness temperature is compared with the ratio of (3,3) to (1,1) 



linewidths in Fig. [13] The ranges of both ratios are divided into 
bins of 0.4 to create the contour plot. We found values between 
0.3 and 5 for the linewidth ratio and from 0.3 to 6 for the (3,3) 
line temperature ratio. The median of Av(3,3)/Av(l,l) is 1.6 with 
a dispersion of 0.8 and the median of Ty[g(3, 3) b s /7MB(3, 3) ca i c 
is 1 .2 with a dispersion of 1.1. Many clumps have higher ob- 
served (3,3) line temperatures than is expected from the calcula- 
tion using the (1,1) and (2,2) line temperatures, which can have 
various reasons: 

1 . Non-thermal excitation can overpopulate or invert the upper 
level of the (3,3) doublet. Different exa mples of these maser 
lines have already been observed: iMangum & Woottenl 
(11994 detected 14 NH 3 (3,3) maser emission towards the 
DR 2 1 (OH) star forming region. iHofner et all ( 1 19941) discov- 
ered the NH3 (5,5) maser line with interferometric observa- 
tions towards the H II region complex G9. 62+0. 19. Some 
ATLASGAL clumps have (3,3) line profiles, which are dif- 
ferent from those of the (1,1) and (2,2) lines and indicate 
maser emission; two examples are shown in Fig. [14J The 
(3,3) line of G25.82-0.18 shows two narrow components 
at lower and higher velocities than v(3, 3) in addition to the 
thermal e mission. This clump also harbours CH3OH and wa- 
ter masers dSzvmczak et al.l2005l) . In the NH3 (3,3) spectrum 
of G30.72-0.08 we see an additional peak at a slightly lower 
velocity than that of the (3,3) thermal emission. This clump 
is located in the W43 high-mass star forming region and is 
revealed as a compact cloud fra gment detected at 1 .3 mm and 
350 //m by iMotte et all (l2003h . However, they did not find 
an association of this source with an OH, H2O or CH3OH 
maser, but instead a bright compact H II region with a flux of 
~ 1 Jy at 3.5 cm. Recent Herschel observations show that it 
is centred on a 1 pc radius IRDC seen in absorpti on at 70 nm 
and b right at wavelengths longer than 160 //m dBallv et al.l 
120101) . However, because we did not find strong maser lines, 
but only weak narrow components in the (3,3) spectra of a 
small subsample of ATLASGAL clumps, it is unlikely that 
population inversion of the (3,3) transition explains the trend 
of increased observed 7"mb(3, 3) values compared to calcu- 
lations. 

2. Since radiative and collisional transitions let the spin orien- 
tations remain constant, transitions between ortho- and para- 
NH3 are not allowed. There are processes a llowing inter- 
conver sion between both species proposed bv lCheung et al.l 
d 19691) . but they are very slow (~ 10 6 yr). The distribution 
betwe en ortho- and para s tates is a strong function of temper- 
ature dTakano et al .120021) : The abundance ratio will be about 
1 if the formation of NH3 occurs in high-temperature regions 
with more than 40 K, but it will rise if NH3 is produced at 
low temperatures, because the lowest level (0,0) is an ortho 
state. Our average r otational temperature is 19 K, for which 
iTakano et alJ d2002l) calculated an ortho-to-para abundance 
ratio of 2. This ratio would result in an increased population 
of the (3,3) level and a decrease of observed to calculated 
(3,3) line temperature ratio. 

3. High rjviB(3, 3) b s values might originate from a hot em- 
bedded source, which is revealed by broad (3,3) linewidths 
leading to a high ratio of Av(3,3) to Av(l,l). We can ex- 
clude any calibration issue, because the NH3 (1,1) to (3,3) 
lines are observed simultaneously and a relative calibration 
error would be the same for all sources, but a wide range 
°f ?Mb(3, 3) b s /7MB(3, 3) ca [ c is found. Since only few typ- 
ical maser line profiles are observed and the ortho-to-para 
ratio would only lead to variations of about a factor two and 
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does not explain the increase in linewidth, we favour the hot 
core explanation. To quantify this effect, we computed the 
resulting rotational and (3,3) temperatures when adding a 
hot core with filling factor i]^ c and temperature of 100 K. In 
Fig. [13] the resulting increase in the (3,3) line temperatures 
is shown for various filling factors. The effect is diminished 
for higher rotational temperatures, since then an observable 
(3,3) line is produced already by the large-scale cold clump 
alone. Sources that are not detected in the NH3 (3,3) line are 
plotted as upper limits (corresponding to three times the rms 
noise in Tmb) m re d- 
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Fig. 12. Width of the NH 3 (3,3) line plotted against that of the 
(1,1) line. The straight line shows equal widths. The contour 
plot is displayed in the top panel with a binning of 0.7 km/s 
for the range of (1,1) linewidths and of 1 km/s for that of (3,3) 
linewidths. 



5. Comparison with the 870 //m dust continuum 

T he ATLASGAL data w ere reduced and calibrated as described 
in lSchuller et alj d2009l) . For the sources of our sample, selected 
as described in Sect. 12.11 we integrated the 870 fim dust con- 
tinuum flux within an aperture with a diameter of 40", which 
corresponds to the beam of the Effelsberg telescope at the am- 
monia inversion frequencies of ~ 24 GHz. The integrated dust 

continuum fluxes S^ 0fjm are given in Table [4] 




Av(3,3)/Av(1 , 




Av(3,3)/Av(1,1) 
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Fig. 13. Ratio of observed to calculated (3,3) main beam bright- 
ness temperatures compared to the ratio of the (3,3) to (1,1) 
linewidths as a contour plot in the top panel and as a scatter plot 
in the middle panel. The correlation plot of the rotational tem- 
perature and observed to calculated (3,3) main beam brightness 
temperature ratio is displayed in the lowest panel. Red points 
show sources that are not detected in the NH3 (3,3) line. Green 
dashed curves indicate very low hot core beam filling factors be- 
tween 0.008% and 0.15%. 



5. 1 . Submillimeter flux and NHj 

The integrated submillimeter dust continuum flux ranges be- 
tween 0.6 Jy and 47 Jy and in this section it is compared with 
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Fig. 14. NH 3 (1,1) to (3,3) lines of two ATLASGAL clumps, 
G25.82-0. 18 at the top and G30.72-0.08 at the bottom: The (3,3) 
line profile is different from those of the (1,1) and (2,2) transi- 
tions, indicating a narrow maser line in addition to the thermal 
emission. 



the ammonia column density. Because the 870 /mi continuum 
emission is optically thin, the submm flux is proportional to the 
dust column density at a given temperature. Assuming a con- 
stant gas-to-dust mass ratio, the dust column density is propor- 
tional to the H2 density (A^h,), which is related to the ammonia 
column density (A^nh,) an d abundance 0fNH 3 ) via the relation 
A^h 2 = A r NH 3 /AfNH 3 ■ Hence, assuming a constant NH3 abundance, 
a correlation between the integrated 870 /mi flux and the ammo- 
nia column density is expected, which is confirmed by this study. 
However, a few clumps exhibit high NH3 column densities and 
low submm fluxes, which can be explained by an increased am- 
monia abundance or an overestimation of the temperature. To in- 
vestigate the two effects we compared the NH3 column density 
with the H2 column density, computed using the NH3 kinetic 
temperatures. 

Models of dust grains with thick ice mantles give an absorption 
coefficient, k, of 1.85 cm 2 /g at 870 /mi at a gas density n(H) = 
10 cm dOssenkopf & Henning|[l99l . We estimated the H2 



column density per the relation dKauffmann et al.ll2008l) 

2.03xl0^-;„/(exp(i^f)-l) 
W(H 2 ) = ^5 V d ; (11) 



k(A)6 2 



Z ' 



where is the 870 /zm flux density in Jy integrated within 

40", A the wavelength in /mi and the dust temperature 
with T(j = 7\i n under the assumption of equal gas and dust 
temperatures. 6 is the 100-m telescope beamwidth of 40" and 
Z/Z the ratio of metallicity to the solar metallicity, assuming 
Z/Z = 1 . We found that most H2 column densities lie between 
2.6 X 10 21 and 1.8 x 10 23 cnr 2 with a peak at 1.7 x 10 22 cnr 2 . 
For each clump we can now use the ammonia column density 
A^nh, to derive the NH3 abundance via^NH^A^NHyA^H,, which 
lies in the range from 7 x 10~ 9 to 1.2 x 10~ 6 with an average 
of ~ 1.2 x 10~ 7 . The H2 column densities, NH3 abundances 
and the integrated dust continuum fluxes X;;„ are given in 

870 j/m 

Table [4] NH3 column densities are plotted against H2 column 
densities in Fig. [15] They show a better correlation than NH3 
column densities and the 870 /mi fluxes, which are affected by 
variations of the temperature of the sources. Since the kinetic 
temperature was used in the calculation of the H2 column 
density, Fig. [15] shows the remaining variation in abundance, 
which can give us a better understanding of the chemical history 
of the clumps. The scatter, given by the ratio of the rms to the 

do" 

mean value, improved from 0.039 for S^ Q m /W N H, to 0.0271 
for NhJNnh^. Deviations from the trend in Fig. [15] are due to 
a higher ammonia abundance, clumps with XNB 3 > 3 x 10~ 7 
are below the lower straight line, or a decreased abundance, 
values lower than 5 xl0~ 8 are above the upper straight line. 
We note that the abundances are upper limits since they are 
rather derived from the ratio of a source and a beam-averaged 
column density, hence affected by the filling factor of the 
cores. With the derived beam filling factors (see Fig. flTT i the 
abundances would on average be higher by about a factor ten. 
With this uncerta inty, we find abunda nces in the range of 5 
x 10~ 9 to 3 x 10~ 7 . lDunham et al.l (1201 ll) found NH 3 abundances 
between 2 x 10~ 9 and 5 x 10~ 7 , wh ich is consistent with our 
results. Moreover, iPillai et all (120061) obtained NH3 abundances 
from 10~ 8 to 10 7 for a sample of IRDCs, their XNH?, range 
is narrower compared to our values and similar to our higher 
abundances. The average of our derived NH3 abundances agrees 
with those predict ed for low-mass pre-pr otostellar cores by 
chemical models dBergin & Langerl 1 1997b . They show that 
ammonia is more abundant in cold dense cores than other 
molecules such as CO, which are depleted from the gas phase at 
high densities. 



5.2. Virial masses and gas masses 

We can estimate the masses of the clumps using parameters de- 
rived from the 870 //m dust continuum and fro m NH3 lines. The 
gas m ass of a source obtained from the dust is dKauffmann et al.l 
12008b 



M gas = 1.2 X 10" 



k(A) 



(12) 



with the distance to the source, d, in pc; for the other param- 
eters see equation [TT] The masses can be determined only for 
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Fig. 15. Correlation plot of the ammonia column density and H2 
column density. Straight lines indicate a high ammonia abun- 
dance ^nh 3 > 3 x 10" 7 , an average value of ~ 1.2 x 10" 7 and a 
low NH3 abundance of 5 x 10" 8 . As the contour plot illustration 
the same correlation is shown in the upper panel, the range of 
NH3 column densities is divided into bins of 10 15 cm" 2 and that 
of H2 column densities into bins of 10 22 cm" 2 . 



the clumps for which we have a preliminary distance estima- 
tion: 277 clumps are associ ated with clouds from the GRS sur- 
vey with known distances dRoman-Duval et a i] |2009b . and 71 
clumps, which are located at the tangential points, have simi- 
lar near and far distances, mark ed in Ta ble |H T o determine dis- 
tances of the clumps studied by iRoman-Duval et aT I (l2009h . we 
followed their choice of near/far distance and used the ammonia 
ve locity together wit h the rotation curve of the Milky Way given 
bv lReid etalJd2009h . The mass can also be calculated using the 
NH3 (1,1) linewidth und er the assumption of virial equilibrium 
dRohlfs & Wilso n 2004) 



M vir = 250Av(l,l) 2 xfl, 



(13) 



with the linewidth Av(l, l) CO n> corrected for the resolution of 
the spectrometer Av(l, l) CO rr = >/Av(L I) 2 _ (0.7km s" 1 ) 2 - The 
radius R in pc is obtained from Gaussian fi ts of the 870 nm 
continuum flux using the Miriad task sfind dSault et al.lll995l) 
at the distance of the source. Equation [13] is derived from the 
virial theorem by neglecting the magnetic energy and the virial 
mass is equal to the dynamical mass. The gas and dynamical 
masses are compared in Fig. [16] the black points in the scatter 
plot indicate sources associated with GRS clouds and the red 



points clumps at the tangential points. Gas masses range from 
60 to 1.6 xlO 4 M Q and virial masses from 20 to 5.5 xlO 3 M Q . 
The straight solid black line shows equal gas and dynamical 
masses. Most data lie below that relation, indicating that for most 
sources their masses exceed their dynamical masses. Since mag- 
netic fi elds can contr ibute to prevent the collapse of molecular 
clouds (ISchulzll2005l) . the magnetic energy has to be taken into 
account in the virial theorem. We did not measure the magnetic 
field and assumed equipartition between magnetic and kinetic 
energy. Then, virialisation is obtained for Mg as =M v i r =2Md vn . 
This relation is shown by the dashed green line, but most sources 
are still located below it . We calculated the virial parameter 
(iBertoldi&McKeell 19921) 



a = 



M, 



(14) 



gas 



which is ~ 1 for clumps supported against gravitational collapse. 
When we use M v ; r = Md yn from equation [13] sources with nar- 
row linewidths have virial parameters on average much smaller 
(mean a = 0.21) than broad linewidth sources with a mean, a, 
of 0.45. With higher density tracers such as H 13 CO + and C 34 S 
on average broader linewidths are found, leading to higher virial 
mass estimates. This will be discussed in Wienen et al. (in prep.). 
Including again the magnetic energy into the virial theorem gives 
a factor 2 resulting in a ~ 1 for clumps with broad linewidth. 
Therefore, these sources are in virial equilibrium. In contrast, 
sources with small Av(l,l) still have a low virial parameter of 
0.42. 

iHill et all d2010h . who observed NH3 lines of high-mass star 
forming regions with the Parkes telescope, found that their sam- 
ple is mostly in virial equilibrium. However, they included many 
sources that have associated infrared sources, UCHII regions, 
and high-mass protostars, but only a few cores in an early evo- 
lutionar y phase w it hout a ny hint at ongoing massive star for- 
mation. IHill et al.1 d2010h found a higher average NH3 (1,1) 
linewidth of 2.9 km s" 1 in contrast to our mean Av(l,l) of 2 
km s" 1 , which results in larger virial masses than we found. 
iDunham et all (1201 Oh obtained an average virial parameter of 
~ 1 from their NH3 observations of high-mass star forming 
sources, which form a similar sample as ours, only at a near dis- 
tance of about 2 kpc compare d to a typical near dist ance of our 
sources of ~ 4 kpc. However. IDunham et al.l d2010l) calculated 
gas and virial masses between 20 and about 1000 M and there- 
fore probed much smaller masses compared to the ATLASGAL 
clumps, whose gas m asses range up to 1.6xl0 4 M . 
IDunham et al.l ( 201 ll) found virial parameters about a factor 2 
larger than our study. Whil e the linewidths and d ust masses in 
both studies are comparable. lDunham et al.l (1201 ll) used for their 
virial mass computation an "effective radius" that on average is 
about twice as large as the half maximum radii that we used from 
ATLASGAL, which accounts for this difference. 
We estimated the magnetic field strength, which is necessary 
to prevent gravitati onal collapse, using B = 2.5 x A^(H2) fiG 
dMcKee et al.ll 1993 1), where NQI2) is the hydrogen column den- 



sity in units of 10 21 cm z . For our values of the hydrogen column 



,-2 



density between 2.3 x 10 and 1.8 x 10 cm"", we obtain mag- 
netic fields between ~ 5 and 450 fiG. 

Estimates of the gas mass and gas column density as calculated 
in Eq. ITTIca n vary strongly. The large uncertainty of Eq . IT2l is 
discussed in lBeuther et al.l d2002l) and lMotte et alj d2007l) , who 
gave a factor 2 of uncertainty in t he mass estimates du e to uncer- 
tain dust emissivity. In addition, iMartin et al.l d2012l) discussed 
evidence for changes of the dust opacity in different environ- 
ments and also found observationally a factor 2 of variation but 
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excluding high-density regions. For denser r egions, direct obser- 
vations of dust opacity changes are difficult. lOrmel et al.l (1201 ll) 
presented new calculations of dust opacities in dense regions as 
a function of grain growth and found for moderate time scales 
(< 10 7 ) a factor 2 increase in the dust opacity with time. The 
effects of the systematic errors of the linewidths that enter Figs. 
[TBI and [16] increase with linewidths, but are much smaller than 
the uncertainties in the dust opacity. As discussed already in 
this section, the main uncertainty here is how representative the 
measured ammonia linewidths are for the dense material of the 
dust clumps. As we will discuss in more detail in Wienen et al. 
(2012b, in prep.), observations of higher density probes towards 
ATLASGAL sources show on average broader linewidths, which 
in turn lead to higher virial mass estimates. 




10 2 10 3 10 4 



M gas/ M 

Fig. 16. Dynamical masses plotted against gas masses for our 
subsamples, which are associated with GRS clouds (black 
points) or located at tangential points (red points). The solid 
black line denotes equal gas and dynamical masses, while the 
dashed green fit shows M gas = 2M<j yn . Note that the gas masses 
are uncertain by a factor 2 due to uncertain dust emissivity. 



6. Correlation of NH 3 and ' CO (1-0) emission 

While ammonia is a high-density tracer (~ 10 4 cm 4 ; 
lUngerechts et alj |1986[) . the 13 CO m olecule can be used t o 
probe low-density regions (~ 10 3 cm~ 3 : ITjngerechts et al.ll 1 997h . 
which are surrounding the dense cores. In this section, the 
ammonia observations are compared with the 13 CO J = 1 — > 
line emission measure at 110.2 GHz in the course of the 
Boston University-Five College Radio Astronomy O bservatory 
Galactic Ring Survey (GRS) dJackson et alj 120061) . It has a 
greater sensitivity of ~ 0.13 K, a higher spectral resolution of 
0.2 km s _1 , comparable or better angular resolution (46") and 
sampling (22") than previous molecular line surveys of the 
inner Galaxy. 13 CO is used instead of the commonly observed 
12 CO J = 1 — > emission, because 13 CO is much less abundant, 
and consequently possesses an optically thinner transition with 
narrower linewidths and its analysis can determine the proper- 
ties of the clumps more reliably. An area of 75.4 deg 2 is covered 
by the whole survey with Galactic longitudes of 1 = 18° - 55.7°, 
limited data are available for 14° < / < 18°, and latitudes of 
\b\ < 1°. The velocity coverage ranges from -5 to 135 km s 
up to a Galactic longitude of 40° and from -5 to 85 km s _1 
for 1 > 40°. We extracted 13 CO (1-0) lines from GRS intensity 



image s, which are integrated over the velocity (Jackson et al. 
120061) . The CLASS software was used to derive the 13 CO LSR 
velocity vo C q of sources detected by GRS and ATLASGAL, 
their linewidths Avi3 C0 and integrated temperatures T MB i3 C0 , 
each with errors. The results are given in Table [5] NH3 and 
13 CO line parameters are compared in correlation plots. 



6.1. Comparison of ammonia and 13 CO lines 

We extracted the 13 CO (1-0) line from GRS intensity images at 
the position of the ammonia observations, fitted a Gaussian to 
the 13 CO line profile and related the 13 CO line with the velocity 
closest to the NH3 velocity to the observed ATLASGAL source. 
Within the ammonia sample we found 517 ATLASGAL clumps 
associated with 13 CO emission. The spectra of 13 CO and NH3 
(1,1) lines are overlaid, some examples of strong sources are 
shown in Fig. [17] The 13 CO linewidths are similar to the widths 
of the NH3 (1,1) inversion transitions for a few sources such as 
G16.91-0.08,G21.39-0.25,G28. 82+0.36, G29.60-0.62, while 
most clumps possess broader 13 CO than NH3 (1,1) linewidths, 
e.g. G15.43+0.19, G18. 82-0.47, G18.88-0.51, G22.37+0.38, 
G18. 83-0.48, and G28.92-0.23. While most ATLASGAL 
clump lines of sight show only one NH3 (1,1) velocity com- 
ponent, the 13 CO (1-0) often shows multiple components (see 
e.g. G16.91-0.08, G21.39-0.25, G28.82+0.36, G15.43+0.19, 
and G22.37+0.38 in Fig.fTTli. Concerning those clumps, we as- 
sumed the 13 CO line with the velocity closest to the NH3 veloc- 
ity to be associated with the ATLASGAL s ource, which has b een 
demonstrated to be a reasonable choice bv lKirk et al.l (120071) for 
low-mass dense cores. 



6.2. Clump-to-cloud motion 

While it is not yet understood how high-mass stars form, two 
main theoretical models are currently under consideration: (1) 
in the turbulent core model massive star forming clouds are sup- 
ported a gainst gravitational collapse mainly by supersonic tur- 
bulence (McKee & "Tknl I2003L iMac Low & Klessenl 120041) . (2) 
Acc ording to recent simulatio ns of star formation in clusters 
(e.g. lBonnell et al.ll997ll200l1) . the accretion rate and in turn the 
mass of the forming star are strongly affected by the density of 
the molecular cloud through which it is moving. In addition, ma- 
terial is supposed to diffuse into the potential well of the cloud at 
the centre of a cluster, where high-mass stars are preferentially 
forming. This theory of competitive accretion also agrees with 
observations, which show th at the most massive sta rs of a clus- 
ter are located at its centre dBontemps et al.ll2010t) . According 
to IWalsh et"aT1 (120041) . this moving accretion is only possible if 
the protostellar core is moving at appreciable velocity with re- 
spect to its surrounding envelope, which will be associated with 
a significant difference in the line-centre velocity of high- and 
low-density tracers. Previous studies of low-mass star forming 
regions have analysed the motions of dense cores with respect to 
their surrounding envelope, using N2FL as high-density tracer 
and 13 CO as wel l as C 1 8 Q as probes of less den se gas (e.g. 
IWalsh etal.ll200l 12007b iKirk et al.l l2007h . While IWalsh et al.l 
(12004 investigated mostly single isolated cores, clustered re- 
gions of dense cores such as NGC1333 and the P erse us molecu- 
lar cloud were also studied bv lWalsh et all (120071) and lKirk et ail 
(120071) . However, they did not find significant motions of the 
cores relative to their lower density surroundings. In this section 
we will compare their results to conditions prevailing in high- 
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Fig. 17. NH3(1,1) inversion transition in red overlaid on the CO(1-0) line in black of a few sources. Some clumps exhibit similar 
widths of the main line of the NH3 (1,1) inversion transitions and I3 CO lines, while most 13 CO linewidths are broader than those of 
the NH 3 (1,1) transitions. 



mass star forming regions. 

Although in the more distant high-mass star forming regions we 
sample clumps on larger scales than the aforementioned stud- 
ies, it is still interesting to investigate motions of clouds rela- 
tive to their host clouds. It is especially useful that the NH3 and 
13 CO lines, good tracers for the former and the latter, respec- 
tively, were measured with similar beamwidths, 40" and 46", 
respectively. The absolute value of the difference in 13 CO and 
NH3 (1,1) line-centre velocities is plotted against the linewidth 
of the ammonia (1,1) inversion transition in the top panel of Fig. 
[T8l If ammonia clumps showed significant motions with respect 
to their envelopes, the differences in line-cent re velocities shoul d 
be comparable to the width of the 13 CO lines (I Walsh et al.l2004l) . 
Deviations in the range of NH3 linewidths would be associated 
with random motions within the ammonia clumps. The dashed 
red line in the top panel of Fig. [18] at 2 km s _1 indicates the av- 
erage widths of NH3 (1,1) lines and the solid green line at ~ 4.5 
km s illustrates the average 13 CO (1-0) linewidths. The rms of 
the difference in both velocities is 0.4 km s _1 , which is shown 
by the dashed dotted black line. Because the distribution of al- 
most all clumps is narrower than the average NH3 linewidth and 
also narrower than the average width of the 13 CO (1-0) transi- 
tion, random motions of the ammonia clumps are expected and 
hint at turbulence, which could prevail in the sources. Hence, 
these low relative velocities between the ammonia clump and 
13 CO cloud put strong constraints on the role of turbulence in 
the clump formation. There are only two sources whose abso- 
lute value of NH3 (1,1) and 13 CO velocity difference is slightly 
higher than the NH3 linewidth: G30.68-0.03 with a deviation in 
both velocities of -2.4 km s has a lower 13 CO velocity, while 
G37.76-0.22 exhibits a lower NH 3 velocity with a difference 



to 13 CO of 2.24 km s" 1 . Fig. [19] illustrates that the 13 CO lines 
of both clumps might show either self-absorption or each 13 CO 
spectrum could be emitted by two clumps, that lie on the same 
line of sight to the observer, leading to higher uncertainties in the 
13 CO velocities. Concerning low-mass star formation, the line- 
centre veloci ty differences of N 2 H + an d 13 CO or C 18 in cores 
observed by IWalsh et al.l (I2004L 120071) are also similar to their 
N2I-F linewidths, which is consistent with the analysis of the 
ATLASGAL sources. This result indicates that the velocity de- 
viations of both molecular probes in low-mass sources are much 
lower than in massive ones, because low-mass cores have a much 
more narrow averag e N2H + and 13 CO l inewidth of ~ 0.38 km s _I 
and ~ 1.18 km s -1 (Walsh et al. 2007) compared to the broader 
ammonia and 13 CO lines of high-mass cores. This also agrees 
with a nd is additionally supported by the analysis of iKirk et al.l 
<l2007h . who measured differences in N->H + and C ls O velocities 
similar to the N2I-F thermal linewidths of 0.065 km s _1 for most 
dense cores in the Perseus molecular cloud assuming a tempera- 
ture of 15 K. 

The distribution of differences in line-centre velocities of 13 CO 
and NH3 in the bottom panel of Fig. [18] illustrates that both, 
higher and lower 13 CO velocities than ammonia velocities are 
present. The differences can be compared to the sound speed of 
the ambient medium, which is 0.23 km s _1 assuming gas that 
consists of 90% H2 and 10% He with a mean molecular weight 
of 2.33 at a temperature of 15 K. Most clumps show differences 
larger than the sound speed, only 32% of them have subsonic 
motions. A high percentage of sources shows a deviation in 
clump-to-cloud velocities between 1 km s and -0.7 km s , 
some range to 2.2 km s and -2.9 km s . This result can be 
compared to the difference in N2H + and C 18 line-centre veloc- 
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ities of low-mass, starless a nd protostellar co res observed in the 
Perseus molecular cloud bv lKirk et alJ (120071) . Although the mo- 
tion of protostellar cores with respect to the envelope could be 
influenced by turbulence generated by heating or outflows, most 
of them as well as the majority of the starless cores show relative 
velocities of both molecules slower than the sound speed, while 
the remaining sources have only slightly highe r differences. This 
agrees with the analysis of IWalsh et al.l ( 120041) . who measured a 
core-to-envelope velocity exceeding the sound speed in only 3% 
of their sample, which included mostly isolated low-mass cores. 
Hence, those comparisons also illustrate that the ATLASGAL 
clumps as a high-mass star forming region sample show much 
more distinct velocity differences between low- and high-density 
tracers than low-mass cores. 
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Fig. 18. Absolute value of the difference in line-centre velocities 
of 13 CO ( 1 -0) and NH 3 lines detected in the ATLASGAL sources 
plotted against the ammonia linewidth in the top panel. The dif- 
ferences of most clumps are similar to the typical NH3 linewidth, 
which is indicated by the dashed red line, and smaller than the 
width of the 13 CO transition labelled by the solid green line, the 
rms of the absolute value of the difference in 13 CO and NH3 
(1,1) velocities is shown by the dashed-dotted line. The bottom 
panel shows the relative number distribution of the ATLASGAL 
clumps with the velocity differences of both molecules. The 
13 CO and NH3 velocities of most sources differ by more than 
the sound speed at 0.23 km s _1 . 
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Fig. 19. NH 3 (1,1) and 13 CO(1-0) lines of G30.68-0.03 are plot- 
ted in the upper panel, the lower one shows these transitions 
of the clump G37.76-0.22. The Gaussian fit of the spectra is 
shown. 



6.3. CO and ammonia linewidths 

The width of the 13 CO (1-0) transition is plotted against that of 
the NH3 (1,1) inversion line in the lower panel of Fig. [20] the 
straight line represents equal widths. For the contour plot in the 
upper panel of Fig. [20] a binning of 0.4 km/s is chosen for the 
NH 3 (1,1) linewidth and of 1 km/s for the 13 CO linewidth. The 
I3 CO line, exhibiting widths from 1 km s _1 to 10 km s _1 with 
a peak at 3.3 km s _1 , is mostly broader than the ammonia line, 
which may indicate that the lines arise from different parts of the 
clumps. It is expected that 13 CO traces a larger scale than am- 
monia, although the beamwidth used for the GRS of 46" is only 
slightly broaderthan that of the Effelsberg telescope. I3 CO traces 
densities of ~ 10 3 cm -3 , but is frozen out at densities of 10 5 
cm -3 , where NH3 is still existing in the gas phase. Hence, am- 
monia probes the massive, denser inner parts of a clump, while 
I 3 CO is also sensitive to the outer less dense regions of the cloud. 
Turbulence within large structures can cause broader 13 CO than 
NH3 linewidths found in small-scale regions of the clumps. 
Another explanation of the difference in linewidths could be that 
multiple I3 CO emitting clumps may lie on the same line of sight 
to the observe r, since 13 CO trac es low densities and its abun- 
dance of 10~ 6 dDu & Yan g 2008) is higher compared to the am- 
monia molecule with 7 x 10~ 9 to 1.2 x 10~ 6 (cf. Sect.EQ. Fig. 
[17]also shows that some ATLASGAL source lines of sight have 
more than one 13 CO component. Those different clumps may 



16 



M. Wienen et al.: Ammonia in massive clumps discovered with ATLASGAL 



also have narrow linewidths, which would add up to the ob- 
served increased 13 CO widths due to velocity dispersion within 
the beam. 




&v(1,1) NH , (km/s) 




Av(1,1) NHj (km/s) 

Fig. 20. Widths of the 13 CO (1-0) transition compared to those 
of the NH3 (1,1) line as contour plot in the upper panel and as 
scatter plot in the lower panel. The straight line corresponds to 
equal widths. The binning used for the NH3 (1,1) linewidth is 
0.4 km/s and for the 13 CO linewidth is 1 km/s. 



7. Temperature and linewidth distributions in 
different evolutionary stages 

To investigate trends of NH3 line parameters with evolutionary 
phases, we divided most clumps into subsamples by searching 
for their associations with other samples that cover the main 
evolutionary stages of high-mass star formation. We used search 
radii dependent on the type of sample. Search radius, number of 
associated sources, mean and error of the NH3 (1,1) linewidth, 
and the rotational temperature of the resulting ATLASGAL sub- 
samples are given in Table [6] While a few sources do not match 
any sample, some clumps show characteristics of more than one 
phase of high-mass star formation and are therefore associated 
with different subsamples. Similar studies using smaller sam- 
ples have already analysed the association of sources traced by 
NH3 with different evolut i onary stages of m assive star forma- 
tion (ISridharan et alJl2005blPillai et alJl2006h . Fig.l2"Tlshows the 



number distribution with the (1,1) linewidth and the rotational 
temperature of the different ATLASGAL subsamples for clarity 
as curves and not histograms, although the data sampling is not 
as good as suggested. 

The lowest panel compares the whole ATLASGAL sample as 
a solid curve in green with a subsample of ATLASGAL clumps 
not detected at 24 //m, denoted as "24 /urn dark" as a dotted curve 
in bla ck, IRDCs as a dashed red curve from iPeretto & Fuller! 
d2009h and a subsample of those IRDCs emitting at 24 /im as a 
dashed-dotted curve in blue. As search radius for associations of 
IRDCs with ATLASGAL source s we used the equivalent radius 
given bv IPeretto & Fuller! d2009l) . which agrees with the radius 
of a disk, whose area is the same as that of the IRDC. An aver- 
age of the search radii used is given in Table [6] The 24 fim dark 
ATLASGAL clumps are the coldest sources with the narrowest 
linewidths, suggesting that they are in a very early phase with- 
out any hint of star formation. The distribution of the IRDCs 
and the ATLASGAL sample is similar in linewidth and rota- 
tional temperature. Although the IRDCs that emit at 24 fim show 
signs of star formation activity, they have only sligthly broader 
linewidths than the average of ATLASGAL sources and simi- 
lar rotational temperatures. The 24 //m emitting star embedded 
in those IRDCs is a local heating component of little influence 
over the volume of the whole clump and thus has no significant 
influence on the source properties. Examples of the NH3 (1,1) to 
(3,3) lines of IRDCs (G18.09-0.30, G10.67-0.22) are displayed 
in Fig. [2] 

We plot the distribution of young ste llar objects (YSOs as dotted 
black curve) from the RMS surve y dHoare et al.ll2005l). CH 3OH 
masers fromlPandian et alj d201 ll) . ISzvmczak et al. d2002l) . and 
ICaswell et al.l (120101) (dashed and dashed-dotted curves) and 
the ATLASGAL sample as a solid curve in green in the mid - 
dle panel. The YSOs are MSX sources dPrice et alj 120011) . 
which emit at infrared wavelengths and are therefore in a more 
evolved phase than the IRDCs. We searched for them within 
the ATLASGAL beam in addition to the MSX positional un- 
certainty. The embedded infrared sources have broader (1,1) 
linewidths than the mean value of ATLASGAL clumps. This 
hints at star formation activity in the YSOs, which results in 
higher rotational temperatures than found for the ATLASGAL 
sample on average (Fig.l2"TTi and at turbulence, e.g. due to molec- 
ular outflows. Concerning the methanol masers, we used two 
search radii for each of the three samples: Sources found within a 
small radius, consisting of the pointing error of the ATLASGAL 
survey of 5" and those of the telescopes used to observe the 
CH3OH samples, might be directly related to the maser (CH3OH 
maser small as a dashed red curve). ATLASGAL clumps lo- 
cated in the same complex as the methanol maser source can 
be revealed within the large radius that includes the Effelsberg 
beamwidth in addition to the pointing errors (CH3OH maser 
large as dashed-dotted blue curve). The range of small and large 
search radii is given in Table [6] The CH3OH samples found 
within both search radii have a similar range of Av(l,l) and 
T I0 u although they have broader linewidths and higher T m i than 
the ATLASGAL sample on average. An ammonia spectrum of a 
YSO (G27.02+0.20) is shown in Fig.|2] 

The upper plot compares the ATLAS GAL sample (the solid 
green curve) with UCHII regions from lHoare et a fl d2005h . We 
again used a small and large search radius as described above to 
distinguish between the embedded object as a dashed red curve, 
and the envelope of the H II region (the dotted black curve). 
Because they already contain a protostar heating up the molecu- 
lar cloud, they are at a more evolved stage of high-mass star for- 
mation than the IRDCs and YSOs and are therefore the warmest 
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sources with the broadest linewidths. Since the embedded object 
traces the central part of the molecular cloud, it covers slightly 
higher temperatures and broader linewidths than the envelope. 
In contrast, the envelope includes cold regions as well, which 
results in a decrease of r ro t and Av(l,l) by averaging over the 
large search radius. Fig.[2]illustrates an example of NH3 lines of 
an embedded UCHIIR (G12.89+0.49). 

To investigate if the differences in physical properties given 
in Table [6] are statistically significant, we performed a 
Kolmogorov-Smirnov (KS) test using the three largest subsam- 
ples from different evolutionary phases. The comparison of the 
NH3 (1,1) line width and rotational tem perature distributions of 
the IRDCs from lPeretto & Fullerl (120091) and the CH3OH sample 
found within the large search radius using the KS test showed 
that both are measurably different (probability > 99.9%). In ad- 
dition, the KS test used for the (1,1) linewidths and rotational 
temperatu res of the IRDCs w ith those of the H II regions enve- 
lope from lHoare et alj d2005l) reveals that those are significantly 
different as well. Cumulative distribution plots (Fig. l22l show 
the NH3 (1,1) linewidth and rotational temperature distribution 
for the different subsamples. They illustrate a trend of broader 
linewidths and higher temperatures for sources with CH3OH 
masers (dashed curve in blue) and with H II regions (dotted black 
line) compared to IRDCs (solid red curve). The distributions of 
methanol masers and H II regions are similar. Our KS test results 
for the NH3 (1,1) linewidths and rotational temperatures do not 
contradict that the CH3OH sample and the H II regions envelope 
sample are drawn from the same parent population. 



8. Conclusions 

The (1,1) to (3,3) ammonia inversion transitions of 898 
ATLASGAL sources, located between 5° and 60° in Galactic 
longitude and ±1.5° in Galactic latitude, were observed with the 
Effelsberg telescope. They build the largest available ammonia 
sample of high-mass star forming clumps. Our results are 
summarised in this section. 

We analysed the correlation of different ammonia line parame- 
ters and compared them to the submm dust continuum and 13 CO 
emission. 

The NH3 (1,1) line's LSR velocity of most clumps lies between 
-10 km s and 120 km s -1 , most clumps sh ow a strong 
correlation with CO emission dDame et alJl200l|) . probing the 
larger giant molecular clouds. The velocities yield near/far 
kinematical distances, that were calculated for 750 unknown 
sources using the revised ro t ation parameters of the Milky 
Way presented by iReid et afl (120091) . Galactocentric distances 
show an enhancement of ATLASGAL sources at 4.5 and 6 
kpc, tracing the Scutum arm and Sagittarius arm, which agrees 
with previous results for o ther tracers of massive star formation 
dAnderson & Baniall2009l) . 

Most sources exhibit NH3 (1,1) linewidths between 0.7 km s _1 
and 3.5 km s _1 , some range up to 7 km s _1 , which is much 
broader than the thermal NH3 l inewidth and com mon linewidths 
observed in cores of low-mass dJiiina etalJl999h . The rotational 
temperature between the (1,1) and (2,2) inversion levels was 
derived and was used to calculate the kinetic temperature. Most 
ATLASGAL clumps have rotational temperatures between 10 K 
and 25 K, some exhibit up to 30 K, while their average kinetic 
temperatures range from 12 K to 35 K, only few up to 40 K. 
The increasing rotational temperature and the broader widths 
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Fig. 21. Relative number distribution of ATLASGAL subsam- 
ples in different evolutionary phases of high-mass star formation 
showing the NH3 (1,1) linewidth in the upper part and the ro- 
tational temperature in the lower part. The whole ATLASGAL 
sample is shown as solid green curve in each panel, while 
the other subsamples are indicated by different patterns and 
colours: Lower panel: dotted black curve: AT LASGAL clumps 
not de tected at 24 /mi, dashed red: IRDCs from lPeretto & Fullerl 
(120091) . dashed-dotted blue: an IRDC subsample emitting at 
24 pm. Middle panel: dotted blac k curve: young st ellar ob- 
jects (YSOs) from the R MS survey dHoare et al.ll2005l). dashed 
red: C H3O H masers fromlPandian et alj ( 201 1 ). Szvmczak et al.l 
(120021) . and lCaswelletal.ld2010l) found within a small search ra- 
dius (see Sect. |7]i, dashed-dotted blue: same CH3OH maser sam- 
ples found in a large radius (see Sect. [7j. Upper panel: UCHII 
regions from iHoare et al.l d2005h . dotted black curve: associa- 
tions with the envelope in a large search radius (see Sect. |7), 
dashed red curve: associations with embedded object in a small 
search radius (see Sect.|7]i. 



of the (3,3) line, which probes higher temperature, hint at star 
formation activity inside the clumps. Indeed, many clumps 
show a stronger than expected intensity of the NH3 (3,3) line as 
estimated for a homogeneous, cold clump with a temperature 
derived from the (1,1) and (2,2) lines. This results probably 
from an embedded hot core component with a low beam filling 
factor and high temperature. 

We used the kinetic temperature and submm dust continuum 
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T„,(1 ,2) (K) 

Fig. 22. Cumulative distribution functions illustrating the NH3 
(1,1) linewidth (upper panel) and rotational temperature (lower 
panel) for the three largest subsamples in T able |6] The distribu- 
tion of IRDCs from lPeretto & Fuller! d2009h is shown as a solid 
red line, the CH3OH sample found within the large search radius 
is plotted as a dashed blue curve and the H II regions envelope 
as a dotted black curve. 

flux to determine the H2 density. Its correlation with the NH3 
column density leads to ammonia abundances relative to H2 
mostly from 5xl(T 8 to 3xlCT 7 . 

For ATLASGAL clumps located at tangential points and asso- 
ciated with clouds from the GRS survey, preliminary distance 
estimates were given and gas masses were calculated, that lie 
in the range from 60 to 10 4 M . Comparing them with the 
virial masses yields a virial parameter of ~ 1 for sources with 
broad linewidths, which are supported against gravitational 
collapse. In contrast, clumps with narrow linewidths have a 
virial parameter of only ~ 0.5. 

The comparison of the NH3 lines as high-density probe with 
the GRS 13 CO emission as low-density envelope tracer yields 
broader 13 CO than ammonia linewidths. This can result from 
turbulence within large structures probed by 13 CO. We also 
analysed the relative motion of the clump and its cloud using 
NH3, which represents the dense core material and 13 CO to trace 
the more diffuse molecular cloud gas. Since our comparison 
yields small differences between NH3 and 13 CO velocities, less 
than the average 13 CO linewidth, random motions of the ammo- 
nia clumps are expected and hint at turbulence. In contrast to 
low-mas s cores, which show mainly subsonic core -to-envelope 
motions dWalshet al.ll2004l 120071: iKirk et al.ll2007l) . the relative 
motions of our sample mostly exceed the sound speed. 
Recent NH 3 (1,1) to (3,3) l ine observations of n orthern BGPS 
sources were carried out by iDunham et al.l d201 lb at the Robert 
F. Byrd Green Bank Telescope (GBT) with a lower FWHM of 



31" compared to the Effelsberg beamw idth. Our results agree 
well with those of lDunham et alJ (1201 lh . yielding similar ranges 
of kinetic temperature, NH3 column density, FL column density, 
thus NH3 abundance, and gas masses. 

Association of the ATLASGAL sources with other high-mass 
star forming samples led to subsamples in various evolutionary 
stages, whose NH3 line parameters were compared. Clumps 
with low temperatures, narrow linewidths, and high column 
densities may be IRDCs, which are assumed to be the cold 
precursors to hot cores and are consequently in an early evo- 
lution ary phase of high-mass star formation dRathborne et al.l 
120061) . Moreover, few ATLASGAL sources show higher tem- 
peratures, greater linewidths and lower column densities than 
the IRDCs. Those properties hint at objects in a more evolved 
stage of massive star formation, such as YSOs that contain a 
luminous infrared source or clumps with embedded UCHII 
regions. 

Because these NH3 observations and analyses can also build 
a valuable database for the selection of different phases of 
massive star formation and the determined velocities are 
important especially for the so far unknown ATLASGAL 
sources, we extended the NH3 survey to the south using the 
Parkes telescope. The comparison of first and fourth quadrant 
properties and the distinction between near and far distances 
will be published in a forthcoming paper (Wienen et al. in prep.). 
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Appendix A: Pointing and calibration 

In this section, we describe the pointing and calibration of the 
NH3 data obtained with the Effelsberg telescope. 
Continuum drift scans were taken on nearby pointing sources 
such as NGC7027, 1741-038, G10.62, and NRAO530 every 
hour to correct for pointing offsets. NGC7027 was also used to 
calibrate the temperature scale. The focal point was determined 
by measuring NGC7027 at the beginning of each observational 
period, before and after sunrise and sunset. 
The data taken with the Effelsberg telescope are given in noise 
tube units and must be converted into main beam brightness tem- 
perature during the calibration procedure. Because many sources 
could only be observed at low elevations between 15° and 25° 
because of their southern position, several flux calibrators at dif- 
ferent elevations had to be observed. These are the continuum 
sources, whose flux density was measured during the pointing. 
In addition, the atmospheric opacity has to be included in the 
calibration procedure owing to the low elevations of the sources. 
The primary flux calibration was performed using NGC7027, 
since its flux density is precisely known. The other pointing 
sources can be used as secondary calibrators, which were ob- 
served over the entire elevation range and had to be adapted to 
NGC7027 as relative calibration. 

The sources were calibrated by the total calibration factor 
Tcal.tob which is given by a constant factor T ca \ belonging to 
the primary flux calibration divided by a gain elevation curve for 
each observation day according to 



(a + b ■ El + c ■ El 2 ) ■ exp(-r • amass) 
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The term yl = a + b ■ El + c ■ El 2 with the elevation El describes 
the distortion of the telescope caused by its weight, when it is 
inclined, as well as thermal conditions. y2= g- Tamass with the 
zenith opacity r and the factor airmass, amass = l/sin(£7X 
takes the weakening of the radiation penetrating the atmosphere 
into account. 

During the primary flux calibration, the factor T ca \ was deter- 
mined using continuum scans of N GC7027, which possesses a 
flux density S A * 5.6 Jy at 1.3 cm dOtt et alj|1994h . The main 
beam brightness temperature Tmb of this cali brator for a given 
beam width 9 and wavelength A is given by dRohlfs & Wilson! 
12001 



(7mb/K) = 



(SjJSy) ■ {A I cm) 2 
(0/arcmin) 2 ■ 2.65' 



(A.2) 



With a beamwidth of 9 — 40" and wavelength of A = 1 .25 cm, 
Tmb of NGC7027 is 7.5 K. We used only one calibration factor 
for the NH3 (1,1) to (3,3) inversion transitions, because the vari- 
ation of the noise diode t emperatur 43 is b etween 4 and 6% and 
of the flux of NGC7027 dOtt et alii 199"! is only ~ 0.1%, both 
are therefore smaller than the average calibration uncertainty of 
10%. 

The opacity is given by the water vapour radiometer, although it 
could not be used in May 2008, because it did not work during 
this month. Hence, the opacity of the sources observed in that 
time was obtained by plotting their system temperatures against 
elevation and fitting them. The system temperature is related to 
the elevation according to 



T^sys - Tq 



7atm(l ~ exp(-r • amass)), 



(A3) 



with the system temperature r sys for an atmospheric temperature 
T a tm, the zenith opacity t and the temperature Tq, which contains 
the temperature of the receiver and the ground. G15.66-0.50 
was observed every day with the Effelsberg telescope. Hence, 
the TjyiB variation of its NH3(1,1) line could be used to estimate 
the rms uncertainty of the calibration, which was determined to 
be ~ 10%, which is small considering the low elevations of the 
sources up to 25°. But this does not influence important NH3 
line parameters such as the optical depth or the rotational and 
kinetic temperature of the clumps, because their determination 
does not require calibration. The optical depth was obtained by 
the ratio of the hyperfine components of one inversion transition. 
The rotational and thus kinetic temperature were calculated by 
the intensity ratio of the ammonia (1,1) and (2,2) lines, which 
can be observed simultaneously. However, the data needed to be 
calibrated to estimate e.g. the source-averaged column density. 



see http://www.mpifr-bonn.mpg. de/div/effelsberg/calibration/1.3cmpf 
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Table 1. NH3(1,1) line parameters. Errors are given in parentheses. The full table is available at CDS. 





RA 1 


Dec 1 


T(U) 


v(l,l) 


Av(l,l) 


Tmb(I-I) 


Name 


(J2000) 


(J2000) 




(km s" 1 ) 


(km s" 1 ) 


(K) 



G5. 39+0.19 17 57 12.9150 -24 12 28.150 1.71 (±0.29) 10.58 (±0.03) 1.42 (±0.06) 1.7 (±0.11) 

G5.35+0.10 17 57 29.2160 -24 16 57.670 2.11 (±0.31) 11.21 (±0.03) 1.48 (±0.07) 1.73 (±0.13) 

G5.64+0.24 17 57 34.6840 -23 58 01.950 2.2 (±0.13) 7.94 (±0.01) 1.61 (±0.03) 4.22 (±0.12) 

G5.62-0.08 17 58 44.6080 -24 08 39.020 3.14 (±0.14) -25.74 (±0.02) 2.43 (±0.03) 3.24 (±0.11) 

G6.22-0.05 17 59 55.2530 -23 36 10.150 3.13 (±0.63) 19 (±0.09) 2.78 (±0.15) 0.62 (±0.09) 

G6.19-0.12 18 00 07.4100 -23 39 51.830 1.5 (±0.41) 10.88 (±0.06) 1.99 (±0.14) 0.99 (±0.12) 

G5. 89-0.29 18 00 07.5280 -24 00 29.480 2.09 (±0.23) 10.11 (±0.02) 1.42 (±0.04) 1.88 (±0.1) 

G5. 89-0.32 1800 14.4580 -2401 23.810 0.9 (±0.13) 10.11 (±0.02) 1.74 (±0.05) 2.6 (±0.09) 

G6.25-0.12 18 00 16.1890 -23 37 09.950 1.71 (±0.25) 14.48 (±0.06) 3.31 (±0.11) 1.02 (±0.08) 

G5.83-0.40 18 00 24.8380 -24 06 50.090 0.91 (±0.21) 7.64 (±0.03) 1.67 (±0.07) 1.54 (±0.08) 

G5.89-0.39 18 00 30.3380 -24 04 00.209 0.69 (±0.09) 9.18 (±0.03) 3.96 (±0.06) 4.05 (±0.13) 

G5.90-0.43 18 00 40.3520 -24 04 18.700 1.8 (±0.06) 6.64 (±0.01) 2.53 (±0.03) 3.8 (±0.09) 

G5. 83-0.51 18 00 51.2040 -24 10 23.290 1.99 (±0.12) 16.22 (±0.01) 1.74 (±0.03) 2.8 (±0.08) 

G6.19-0.36 18 01 01.7140 -23 47 05.130 2.11 (±0.12) -33.81 (±0.02) 2.2 (±0.04) 3.72 (±0.12) 

G5.91-0.54 18 01 08.2520 -24 07 14.800 2.09 (±0.22) 14.79 (±0.02) 1.37 (±0.05) 2.32 (±0.11) 

Notes. (1) Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
Table 2. NITi(2,2) and (3,3) line parameters. Errors are given in parentheses. The full table is available at CDS. 





v(2,2) 


Av(2,2) 


T MB(2,2) 


v(3,3) 


Av(3,3) 


T MB (3,3) 


Name 


(km s" 1 ) 


(km s" 1 ) 


(K) 


(km s- 1 ) 


(km s 


') 


(K) 


G5.39+0.19 


10.6 (±0.13) 


1.53 (±0.38) 


0.47 (±0.12) 








< 0.09 


G5.35+0.10 


10.85 (±0.17) 


1.8 (±0.35) 


0.38 (±0.11) 








<0.08 


G5.64+0.24 


7.95 (±0.02) 


1.89 (±0.06) 


2.88 (±0.12) 


8.29 (±0.09) 


2.81 (J 


:0.2) 


0.82 (±0.09) 


G5.62-0.08 


-25.57 (±0.05) 


2.93 (±0.11) 


1.93 (±0.14) 


-25.47 (±0.05) 


3.54 (i 


=0.13) 


1.85 (±0.11) 


G6.22-0.05 






<0.1 








<0.08 


G6.19-0.12 


10.63 (±0.26) 


2.91 (±0.62) 


0.3 (±0.1) 








<0.1 


G5.89-0.29 


10.13 (±0.06) 


1.65 (±0.16) 


0.85 (±0.1) 








<0.1 


G5.89-0.32 


10.09 (±0.06) 


1.87 (±0.14) 


1.21 (±0.13) 


10.04 (±0.06) 


1.51 (d 


:0.23) 


1.01 (±0.1) 


G6.25-0.12 


13.97 (±0.2) 


3.02 (±0.42) 


0.44 (±0.12) 


14.35 (±0.39) 


4.06 (d 


= 1.18) 


0.31 (±0.1) 


G5.83-0.40 


7.77 (±0.11) 


2.48 (±0.3) 


0.67 (±0.11) 








<0.09 


G5.89-0.39 


9.28 (±0.04) 


4.58 (±0.09) 


3.37 (±0.15) 










G5.90-0.43 


6.58 (±0.04) 


3.29 (±0.08) 


2.92 (±0.14) 


6.6 (±0.05) 


3.87 (d 


:0.12) 


1.76 (±0.09) 


G5.83-0.51 


16.16 (±0.06) 


2.35 (±0.14) 


1.13 (±0.1) 
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G6.19-0.36 


-33.84 (±0.04) 
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2.71 (±0.14) 
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1.82 (±0.09) 
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Table 3. Parameters derived from the NH3(1,1) to (3,3) inversion transitions. Errors are given in parentheses. The full table is 
available at CDS. 







Trot 






T kin 


log 10 (N NHi ) 

(CITT 2 ) 


Name 


(K) 






(K) 


G5 


39+0.19 


13.12 (d 


1 


27) 


14.39 (±1.67) 


15.12 (±0.08) 


G5 


35+0.10 


11.66 (d 


1 


15) 


12.52 (±1.43) 


15.23 (±0.07) 


G5 


64+0.24 


19.53 (d 





82) 


23.94 (±1.39) 


15.34 (±0.03) 


G5 


62-0.08 


16.08 (d 





81) 


18.5 (±1.19) 


15.63 (±0.02) 


G6 


22-0.05 


- 






- 


- 


G6 


19-0.12 


13.92 (d 


1 


99) 


15.45 (±2.7) 


15.21 (±0.12) 


G5 


89-0.29 


15.45 0* 


1 


02) 


17.58 (±1.47) 


15.21 (±0.05) 


G5 


89-0.32 


18.16 (d 


1 


08) 


21.67 (±1.73) 


14.97 (±0.07) 


G6 


25-0.12 


15.79 (d 


2 


13) 


18.07 (±3.09) 


15.5 (±0.07) 


G5 


83-0.40 


17.51 (d 


1 


53) 


20.66 (±2.38) 


14.95 (±0.1) 


G5 


89-0.39 


27.33 (d 


1 


33) 


39.65 (±3.21) 


15.36 (±0.06) 


G5 


90-0.43 


22.69 (d 


1 


15) 


29.63 (±2.22) 


15.5 (±0.03) 


G5 


83-0.51 


14.82 (d 





64) 


16.68 (±0.89) 


15.28 (±0.03) 


G6 


19-0.36 


20.85 (d 


1 


17) 


26.22 (±2.08) 


15.48 (±0.03) 


G5 


91-0.54 


17.44 (d 


1 


22) 


20.55 (±1.89) 


15.22 (±0.05) 



Table 4. Parameters from the 870 /urn dust continuum and NH3 lines. The full table is available at CDS. 









X 


o40" 

°87(Vm 


Distance (near) 


Distance (far) 


Name 


1 


(10 22 cm" 2 ) 


(10- 7 ) 


(Jy) 


(kpc) 


(kpc) 


G5.39+0.19 


0.18 


1.18 


1.12 


1.21 


2.8 


13.9 


G5.35+0.10 


0.2 


1.79 


0.95 


1.44 


2.9 


13.8 


G5.64+0.24 


0.22 


3.79 


0.58 


8.42 


2.2 


14.5 


G5.62-0.08 


0.21 


2.59 


1.65 


3.96 






G6.22-0.05 








2.74 


3.7 


13.0 


G6.19-0.12 


0.1 


1.32 


1.23 


1.52 


2.7 


14.1 


G5.89-0.29 


0.14 


2.1 


0.77 


2.98 


2.5 


14.2 


G5.89-0.32 


0.23 


2.15 


0.43 


4.15 


2.6 


14.1 


G6.25-0.12 


0.08 


2.16 


1.46 


3.19 


3.1 


13.6 


G5. 83-0.40 


0.14 


1.45 


0.61 


2.62 


2.1 


14.6 


G5.89-0.39 


0.22 


9.74 


0.24 


41.64 


2.4 


14.4 


G5. 90-0.43 


0.17 


6.88 


0.46 


20.37 


1.9 


14.8 


G5.83-0.51 


0.23 


3.8 


0.5 


4.96 


3.4 


13.3 


G6.19-0.36 


0.18 


3.73 


0.81 


9.38 






G5.91-0.54 


0.15 


1.37 


1.21 


2.45 


3.3 


13.5 



Notes. Columns are beam filling factor, H 2 column density, NH 3 abundance, 870 yum dust continuum flux integrated over 40", near kinematic 
distance, far kinematic distance. * are sources with known distances (cf. Sect. 15. 21 ) 

Table 5. Line parameters derived from 13 CO (1-0) transition. Errors are given in parentheses. The full table is available at CDS. 



Name 


Vl3 CO 
(km s 






Av,3 co 
(kms- 1 ) 


T MB. 
(K) 


3 CO 


G14.62+0.33 


26.73 


(±0 


03) 


2.56 (: 


t0.09) 


3.26 (: 


tO 


17) 


G14.63+0.31 


26.03 


(±0 


04) 


2.34 (: 


tO.l) 


3.65 (: 


tO 


16) 


G14.92+0.07 


25.93 


(±0 


05) 


3.35 (: 


t0.12) 


3.76 (: 


tO 


26) 


G15.01+0.01 


25.38 


(±0 


13) 


4.35 (: 


t0.42) 


2.65 (: 


tO 


15) 


G15.43+0.19 


47.85 


(±0 


03) 


3.13 (: 


t0.06) 


3.31 (: 


tO 


17) 


G14.99-0.12 


49.46 


(±0 


06) 


3.15 (: 


t0.19) 


1.67 (: 


tO 


18) 


G15.46-0.00 


66.92 


(±0 


06) 


3.3 (± 


0.13) 


1.66 (: 


tO 


12) 


G14.89-0.40 


61.64 


(±0 


13) 


2.65 (: 


t0.32) 


2.74 (: 


tO 


66) 


G15.20-0.44 


20.86 


(±0 


03) 


2.38 (: 


t0.08) 


3.69 (: 


tO 


13) 


G15.53-0.41 


39.27 


(±0 


07) 


4.36 (: 


t0.18) 


1.97 (: 


tO 


19) 


G15.56-0.46 


15.24 


(±0 


03) 


1.78 (: 


t0.06) 


3.92 (: 


tO 


28) 


G16.36-0.07 


47.25 


(±0 


05) 


2.72 (: 


t0.13) 


5.75 (: 


tO 


19) 


G16.58-0.05 


59.41 


(±0 


02) 


3.27 (: 


t0.06) 


6.53 (: 


tO 


17) 


G16.83+0.08 


62.43 


(±0 


03) 


2.48 (: 


t0.07) 


2.35 (: 


tO 


15) 


G16.58-0.08 


40.85 


(±0 


17) 


3.74 (: 


t0.57) 


2.32 (: 


tO 


13) 
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6. Properties of ATLASGAL subsamples belonging to different evolutionary stages of high-mass star formation. 



Sample 


Search radius (") 


Number 


Trot (K) 


o-(Trot) 


Av(l,l)(kms _1 ) 


o-(Av(l,l)) 


24yum dark 


40 


98 


15.25 


0.21 


1.85 


0.09 


IRDCs 


32.27 


264 


16.31 


0.19 


1.96 


0.04 


IRDCs24yum 


32.27 


117 


16.89 


0.27 


2.08 


0.06 


YSOs 


40.05 


57 


17.31 


0.35 


2.14 


0.08 


CH3OH maser small 


5-25 


61 


19.23 


1.61 


2.77 


0.1 


CH3OH maser large 


40-47 


137 


18.97 


1.65 


2.65 


0.1 


HII regions embedded object 


5 


57 


18.98 


0.37 


2.79 


0.18 


HII regions envelope 


40.3 


194 


18.63 


0.2 


2.63 


0.08 



